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We present the first calculations of the electromagnetic form factors of 3 He and 3 H within the 
framework of the Covariant Spectator Theory (CST). This first exploratory study concentrates on 
the sensitivity of the form factors to the strength of the scalar meson-nucleon off-shell coupling, 
known from previous studies to have a strong influence on the three-body binding energy. Re- 
sults presented here were obtained using the complete impulse approximation (CIA), which includes 
contributions of relativistic origin that appear as two-body corrections in a non-relativistic frame- 
work, such as "Z-graphs" , but omits other two and three-body currents. We compare our results to 
non-relativistic calculations augmented by relativistic corrections of 0(v/c) 2 . 



I. INTRODUCTION 

The Covariant Spectator Theory (CST) 0, 0, [| is 
a manifestly covariant formalism developed in particu- 
lar for the description of few-nucleon systems. One of 
its most characteristic features is that four-momenta are 
conserved in intermediate states and all particles but one 
are confined to their mass shells. The remaining particle 
is off mass shell with its propagation described by the 
corresponding full Feynman propagator. One can derive 
a closed set of integral equations for the scattering am- 
plitude, which effectively sum an infinite set of Feynman 
diagrams that consist of iterations of a kernel with an off- 
shell propagator after the on-shell constraints have been 
incoporated. The kernel can be symmetrized (or anti- 
symmetrized) as required for the treatment of identical 
particles. 

For the nucleon-nucleon (NN) scattering amplitude, a 
number of CST potential models of One-Boson-Exchangc 
(OBE) form have been developed that give a good de- 
scription of the deuteron and the elastic NN scattering 
observables below 350 MeV 0, H, @. The presence of 
off-mass-shell nucleons allows for a richer structure of 
the OBE kernels than what is possible in nonrelativistic 
frameworks. The general structure of the CST kernels 
contains terms for the coupling of bosons to nucleons 
that are exactly zero when the nucleon is on mass shell 
but contribute for off-mass-shell nucleons. For the scalar- 
isoscalar (a or ctq) and scalar-isovector (6 or a±) mesons 
these couplings were investigated for the first time in Q , 
and it was found that they increase the quality of the fit 
to the observables significantly. 

When the CST equation for the three-nucleon (3N) 
bound state [|| was solved for the first time, it was real- 
ized that the scalar off-shell coupling strongly influences 
the 3N binding energy 0] • When the scalar off-shell cou- 
pling strength was varied systematically, it turned out 



that the value that gave the best fit to the NN ob- 
servables simultaneously also produced the best agree- 
ment with the experimental triton bound state energy of 
E t — —8.48 MeV. What is also remarkable about this 
result is that nonrelativistic calculations with so-called 
"realistic" potentials, i.e., potentials that fit the NN scat- 
tering data with a x 2 /-^data ~ 1, do not reproduce E t 
without the addition of irreducible 3N forces specially 
adjusted for that purpose. 

This somewhat unexpected success of the CST NN po- 
tential models in explaining E t appeared at first to be to 
some extent accidental. However, recently we found two 
new high-precision OBE models that fit the world data 
on np scattering below 350 MeV with an essentially per- 
fect x 2 /A r data ~ 1 [3, @], and in both cases we again 
found the triton binding energy to be very close to the 
measured value. Thus we are lead to conjecture that our 
relativistic OBE kernels capture the essential part of the 
physics responsible for the binding of the 3N system. 

Of course one expects more from a realistic descrip- 
tion of a bound state than merely to reproduce the cor- 
rect binding energy, ft should also reproduce its internal 
structure, which can be accessed through the correspond- 
ing electromagnetic form factors. While the interaction 
of nuclei with electromagnetic probes introduces some 
new problems in the form of ambiguities in the definition 
of the nuclear currents, the calculation of the electromag- 
netic 3N form factors certainly presents an interesting 
and necessary test of the three-body CST. 

The exact form of the electromagnetic currents of a 3N 
system in CST was derived in a previous paper Q. In 
the present work we use this current in order to calcu- 
late the electromagnetic form factors of 3 H and 3 He, as 
well as their isoscalar and isovector combinations, in the 
so-called complete impulse approximation (CIA), i.e., in 
an approximation that includes all one-body currents but 
leaves out two- and three-body currents. As will be dis- 
cussed, our CIA includes, in principle, contributions that 
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in nonrelativistic frameworks appear as interaction cur- 
rents, so care has to be taken when results of different 
formalisms are compared. 

It is well known that interaction currents can have a 
significant influence on the electromagnetic 3N form fac- 
tors. Therefore, a calculation in CIA cannot be expected 
to provide a quantitative description of the data over the 
entire range of the transferred momentum. Neverthe- 
less, there is a number of important aspects that can be 
learned from a CST calculation in CIA: 

(a) The first aspect is technical in nature. The calcula- 
tion of the form factors requires the numerical computa- 
tion of rather envolved expressions, some of which have a 
complicated singularity structure. Since this is the first 
time these form factors are calculated in CST, it is im- 
portant to learn how this can be done in an efficient and 
numerically stable way. In the same category belongs the 
study of the convergence of the employed expansion into 
hclicity partial waves. 

(b) For a given 3N vertex function obtained by solv- 
ing the CST 3N bound state equation with a chosen 
NN interaction model, we study the relative importance 
of the various different contributions to the total re- 
sult. For instance, there are six Feynman diagrams that 
define the CIA, and it is useful to determine whether 
they are of comparable magnitude or if one may neglect 
some of them. Another interesting question is the im- 
portance of negative-energy states. The Feynman prop- 
agator of an off-shell nucleon can be decomposed into 
positive- and negative-energy parts, or equivalently into 
parts with positive or negative p-spin. The latter are re- 
lated to the "pair terms" or "Z-graphs" in the language 
of time-ordered perturbation theory. It is known that, 
in nonrelativistic calculations with relativistic correc- 
tions the Z-graphs connected to one-pion exchange have 
a sizeable effect. While a direct quantitative compar- 
ison of our negative-energy-state contributions to these 
pion Z-graphs is not easy because the decomposition into 
positive- and negative-energy states is frame dependent, 
one may still get a general idea about their size. 

(c) Still operating with one particular 3N vertex func- 
tion, we perform sensitivity studies with respect to the 
parameterization of the single- nucleon current. In addi- 
tion to changes in the Dirac and Pauli nucleon form fac- 
tors, off-shell nucleon currents also have off-shell terms 
and their associated form factors. We investigate the 
sensitivity of the results on changes in these off-shell cur- 
rents. 

(d) We investigate the model dependence of the form 
factors by varying the NN interaction model from which 
the 3N vertex function was computed. This is done by us- 
ing various members of a family of potentials that were 
fitted with different, fixed values of the scalar off-shell 
coupling. These models produce rather different 3N bind- 
ing energies while the \ 2 to the data varies only mod- 
erately, which allows one to assess also binding energy 
effects on the form factors in a systematic way. 

(e) Since the calculations in the present work use NN 



interaction models in which the pion-nucleon coupling is 
purely of pseudo-vector form, the coupling to negative- 
energy nucleon states (or Z-graphs) is expected to be 
suppressed. Therefore it makes sense to compare our re- 
sults to the ones in impulse approximation obtained in 
a nonrelativistic framework [?|, where relativistic effects 
are calculated to order (v/c) 2 , while in CST they are in- 
cluded in all orders. The comparison may thus help to 
decide whether the (v/c) expansions of relativistic cor- 
rections are justified. 

(f) Last, but certainly not least, the calculations in this 
work give an indication of whether or not the CST cal- 
culations of three-body form factors will require any new 
physics. For example, when we first calculated the three- 
nucleon binding energies in the CST we obtained about 
—6 MeV, much too underbound. In trying to under- 
stand this result we discovered that the binding energy 
was very sensitive to the off-shell couplings of the scalar 
mesons, and that the introduction of such couplings not 
only "corrected" the binding energy but also improved 
the fits to the two-body data. We might not have real- 
ized that such previously not considered couplings were 
needed to give a good fit to the two-body data if we had 
not been initially confronted with such a large discrep- 
ancy in the binding energy. In short, solving the binding 
energy "problem" lead to the discovery of a new inter- 
action that, in CST, is also necessary for an efficient de- 
scription of two-body data. In this context, study of the 
form factors not only serves as a stringent test of these 
off-shell scalar couplings, but also exposes us to the pos- 
sibility that another new mechanism will be required for 
the understanding of the form factors. We return to this 
issue at the end of the paper. 

This paper is organized in the following way: after 
the introduction in section I, section II provides a brief 
overview of the formalism in which the electromagnetic 
3N form factors were calculated. The numerical results 
are presented and discussed in section III, and section 
IV contains our conclusions. Appendices A and B show 
some details about how the calculations were carried out. 



II. FORMALISM 

A. Charge and magnetic form factors of composite 
spin-1/2 particles 

We start by reminding the reader that the most general 
form of the spin 1/2 current of the 3N bound states (with 
mass M t ) depends on two form factors 



(M'\J^\M) = u{P' t ,M<) Ff N (Q 2 ) 7 P 



+F* N m 



2M t 



u(P t ,M), (1) 



where q = P[ — P t is the photon momentum, P t and P[ 
arc the initial and final trinucleon momenta, M and M' 
the respective spin projections along the z-axis, , and we 
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use the convention Q 2 = — q 2 . (Note the use of the roman 
q to represent the photon four- momentum, instead of the 
more customary q.) 

It is convenient to rewrite this current using the Gor- 
don decomposition (valid between on-shell spinors), 



2M f 



(Pl + PtY 
1M t 



(2) 



and to evaluate it in the Breit frame, where q° = and 
P' t = q /2 = -P t . We obtain 

(M'\J° m \M) 



S(P t ',M')[F A 3 f(Q 2 ) 7 ° 
E 



-f 2 (Q 



M t 



u{P u M), 



(M'\r 3N \M) = F 3 ?{Q 2 )u{P^M' 



2M t 



(3) 



u{P u M), 



(4) 



where E = P/° = P t ° and F M = F 1 + F 2 is the familiar 
magnetic form factor. Equation (j4]) shows immediately 



that, if we choose q to lie in the +z direction, J 3N = 
and the spatial components of the current depend only 
on the magnetic form factor. 

To evaluate these matrix elements we first start with 
the spinors for the three-body states when they are at 
rest: 



u(0,M) 



where M = ±5 with 



XI = 



X-i = 



(5) 



(6) 



Then, the initial and final three-body states will be ob- 
tained by boosting the rest spinors to the correct frame. 
The Dirac operator for an active boost to a frame with 
velocity sinh£ = Q/2M t in the +z direction (we use 
Q = |q|) is denoted S(B(£,e 3 )) and defined in Eq. @. 

With this notation for the Breit frame matrix elements, 
where the initial state is boosted in the — I direction, the 
matrix element (|3|) becomes 



Breit frame conserves spin (or flips helicity) and equals 
the charge form factor F c = F 1 -tF 2 with t = Q 2 /AM 2 . 
Similarly, the x component of the current is 



{M'\Jl N \M) 



M 

2M t 



■«(0,M')5- 1 (B(fe 3 )) 
QF ^ r ^u(0 ) M')iE 2 u(0,M) 



2M t 



= 2M'S 



M' , —M ' 



Q F lf{Q 2 ) 



M 

2M t 



(8) 



where we replaced — ia 13, = iY, 2 [with the £j matrices in 
Dirac space defined in Eq. (|27p] and the phase 2M' arises 
from the matrix element of iS 2 . 

The y component of the current is obtained in the same 
way, yielding 



(M'\Ji N \M)=i2M'6 M ,,- M - 



Flf{Q 2 ) 



2M t 



(9) 



In this work, we actually calculate the 3N current matrix 
elements in the Lab frame. Since the current matrix ele- 
ments transform like a four-vector, we can apply a boost 
from the Breit to the Lab frame, 

(MVL>\M) = [B-mW (M'\J£ ieit \M) , (10) 

in order to determine the relation between the charge and 
magnetic form factors and the Lab frame current matrix 
elements. 



Fc N {Q 2 
Fl?{Q 2 ) = -(2M) 



\- -21 

2M t , 



i(2M) 



2M, 



-M\Jl N ^ h \M) 



-M\jl N r ab \M). 



(11) 

(12) 
(13) 



We conclude that the charge and magnetic 3N form fac- 
tors can be obtained from the matrix elements of the time 
and space components of the 3N electromagnetic current 
in any convenient frame. 



(M'\J° 3N \M) = f2(0,Af')^ 1 ( J B(ee 3 ))^M(g 2 )7° 
S(B(-te 3 ))u(0,M) 



F 2 {Q 2 )^- 
W ' M t 



= £i(0, M') 

= $M,M 



F M (Q 2 )-F 2 (Q 2 



M 2 



u(0,M) 



Fi(Q 2 )~^F 2 (Q 2 ) 



(7) 



where E Q = M 4 (cosh 2 ^/2 + sinh 2 £/2) = A/ t cosh^ = 
\J M 2 + Q 2 /4. The 0th component of the current in the 



B. The three-nucleon vertex functions and the 
relativistic wave function 

The three-body form factors in the spectator theory 
are expressed in terms of the bound-state Faddeev ver- 
tex function T\ 1 x 2a (ki, k 2 , k 3 ) with (by convention) par- 
ticles 1 and 2 on mass shell with helicities Ai and A2, 
and particle 3 off mass shell with the associated Dirac 
index a {ki is the four-momentum of nuclcon i with 
mass m) . The total four-momentum (which is conserved) 
is P = ki + k 2 + k 3 and k\ — k 2 = to 2 . The on- 
shell energy of a nucleon with momentum ki is denoted 
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FIG. 1: (Color online) Diagrammatic representation of the Fad- 
deev equation dlTI . 



FIG. 2: (Color online) Diagrammatic representation of the 
Eq. i21l for the symmetrized two-body scattering subamplitudc 



with both final particles off shell. 



E(ki) = y/rn 2 + k|. In this paper we adopt the alter- 
native notation for T = T(ki, k 2 \ P), convenient because 
k 3 = P — k\ — k 2 is a dependent variable. This vertex 
function is related to the relativistic wave function ^, 
defined by the relation 

^\ 1 \ 2a {k u k 2 ;P) = G QQ /(fc 3 )r AlA2Q /(fci,/c 2 ;P), (14) 

where we continue to use k 3 as the argument of G to 
make the notation more compact. The propagator G is 

1 



G aa > (k 3 



(to - # 3 
m \ - 



E(k 3 



A 



u a < (k 3 , A)u Q (k 3 , A) 



« 3 



E{k 3 ) - ie 
v a /(-ks, X)v a (-k 3 , A) 



^3 



E(k 3 ) 



E(k 3 



Xp 



/9 <,(k 3 ,A)<(k 3 ,A) 



^3 



pE{k 3 



it 



(15) 



where E(k 3 ) is the on-shell energy of particle 3, fc 3 its 
off-shell energy, and the p-spin is either 4- or — with the 
convention 



w„(k 3 ,A) 
"a(k 3 ,A) 



it a (k 3 , A) 
v a (-k 3 , A) . 



(16) 



The vertex function satisfies the following Faddeev 
equation (shown diagrammatically in Fig.[TJ with £ = +1 
for bosons and C = — 1 for fcrmions) 

r lu u v\ ■ md3k! > 

L\ 1 \ 2 a\k 1 ,k 2 ;P) 



E(k' 2 ) (2tt)3 
,{k 2 ,k' 2 ;P 23 ) 2(V 12 *a 1 aSo'(Ai, k' 2 ; P) , 

(17) 

where P 23 = P — k± is the total four-momentum of the 
pair, M is the two-body scattering amplitude satisfying 

M\ 2a ,\< 2a > (k 2 , k' 2 ; P 23 ) = V\ 2a ,\> 2a >(k 2 ,k' 2 ;P 23 ) 



md 3 k 2 



Y. V WH3{k 2 ,k'i;P 23 ) 



A" 



xG^(P 23 - k 2 ')M x ,,p, >x , ,(k 2 \ k' 2 - P 23 ) , (18) 



V\ 2 is the permutation operator that interchanges par- 
ticles 1 and 2, and V is the two body interaction kernel 
(which we will frequently call the "potential" because of 
its close connection with the nonrelativistic potential). 
For convenience we have adopted the notation 

V\ 20lt x 2 a'(k 2 , k' 2 ;P 23 ) 

= Uf3(k 2 ,X 2 )V0 a ,0' a '(k 2 ,k 2 ;P 23 )up'(k 2 ,X' 2 ) (19) 

so care must be taken to distinguish Dirac indicics from 
helicity indices. Whenever a Dirac index is replaced by 
a helicity index, a contraction with an on-shell, positive 
energy spinor, such is shown in Eq. (|19p . is implied, and 
unless otherwise stated, it is assumed that the particle 
is on-shell. The on-shell Dirac spinors are normalized to 
uu = 1. 

Calculations of the form factors requires knowledge of 
the vertex function with the two interacting nucleons off- 
shell. This vertex function was defined in Ref. [8j and can 
be obtained using the Faddeev Eq. ifTT)) . generalized to 
the case when both of the final state interacting nucleons 
are off-shell 



^Xifiaiki, k 2 ;P) = - J 



E(k' 2 )(2nf 



x M 0oc,x' 2 a'(k 2 , k' 2 ;P 23 )2(V 12 ^ Xl X' 2 a'(ki, k' 2 ; P) , 

(20) 

where now k 2 ^ m 2 . The off-shell scattering amplitude is 
obtained by quadratures from the off-shell kernel Vp at x> a ' 
and the on-shell scattering amplitude M\ iipi \> a > using a 
generalization of the two-body equation (fhSf 

Mp at x 2 a'(k 2 ,k 2 ;P 23 ) = Vp a ^' 2 a'{k 2 ,k' 2 ;P 23 ) 



E(k>>) (2tt)3^ 



xGV(P 23 - k'{)M K p,^ a ,{k' 2 \k' 2 -P 23 ) . (21) 

This equation is illustrated in Fig. [2j The off-shell kernel 
V/3a,X' a' is known (in principal), and is discussed in more 
detail in Appendix [XJ 



C. The three- nucleon form factors 

The diagrams needed to calculate these form factors 
are displayed in Fig. [3j The diagrams A-F are referred to 
as the complete impulse approximation (CIA). Diagrams 
G-J are the interaction currents. The algebraic result for 
the six diagrams that make up the CIA can be written 



X 



(A) ^ 



+ 6? 





+12 



FIG. 3: (Color online) The electromagnetic 3N current in CST for elastic electron scattering from the 3N bound state. Diagrams (A) to 
(F) define the complete impulse approxiumation (CIA), in which the photon couples to single nucleons, which can be off-shell (A and D) 
or on-shell before or after the photon- nucleon vertex (B, C, E, and F). The interaction diagrams (G) to (J) describe processes in which 
gthc photon couples to two-body currents associated with the two-nucleon kernel. 



J CIA 



= 3e 



m 2 d 3 kid 3 k 2 



he E {^AiW (h, ha; Pi) [1 + 2 CPu] %M,te) ^x 1 x 2a (h,k 2 ;P t ) 



Ai A2 



E{kx)E(k 2 ) (2i 

+F Al 0'a(fci, k' 2 ; P>) Gp^{k' 2 )j^{k' 2 ,k 2 ) u 7 (k 2 , A 2 ) [1 + 2 £V 12 ] * AlAa a(fci, k 2 ;P t ) 
+^x 1 x 2a (k 1 ,k' 2 ;Pl)[l + 2 CV 12 ] u 7 (k' 2 , X 2 ) jfc (4 , k 2 ) Gpp (k 2 ) T XlPa (fc x , k 2 ; P t ) } , 

I 



(22) 



where the doubly off-shell vertex functions are evaluated 
using ([20]) . j a 'a(k' , k) is the single nucleon current for off- 
shell nucleons with incoming (outgoing) four- moment a 
fcn {k'^ 1 ), Pt is the four-momentum of the incoming, and 
Pi = P t + q of the outgoing three-body nucleus, k[ = 



ki + q, and in every term kf = k. 



Each diagram 



has two off-shell propagators and three interactions, for 
a phase of i(—i) 5 — 1. Mathematical expressions for 
the diagrams with interaction currents are given in|8| 
(similar expressions are also derived in 

and [llf), 

and will not be needed here. Equation (|2"2"]) with the 
two particle off-shell vertex function defined by PD)) is 
evaluated numerically in this paper. 



D. Three- nucleon partial wave helicity states 

In the 3N bound states, the three nucleons are treated 
as identical particles with an isospin degree of freedom. 
The expression for the 3N current should therefore be 
symmetric with respect to which nucleon the photon cou- 
ples to, and which nucleons are on or off mass shell. How- 
ever, permutation symmetry allows us to relate all pos- 
sible diagrams that contribute to the current to ones in 
which nucleon 1 is a spectator and on mass-shell, while 
nucleons 2 and 3 form a pair that undergoes a two-body 



interaction before or after the photon couples to one of 
them (in CIA). The photon can also couple to the inter- 
acting pair directly, which produces a two-body current. 
The resulting total 3N current is displayed in Fig. [3] 

A closer inspection of this 3N current shows that par- 
ticle 1 is always on mass-shell, particle 3 always off mass- 
shell, and particle 2 in some cases on and in others off 
mass-shell. It is therefore sufficient to introduce 3N basis 
states with these characteristics when we want to express 
amplitudes as matrix elements of operators between 3N 
states. 

The propagator of an off-shell particle can be decom- 
posed into positive- and negative-energy contributions, 
as described by Eq. (fT5]) . where also the p-spin notation 
was introduced. In terms of the basis states, the presence 
of an off-shell particle reflects itself in the appearance of 
two p-spin states p = ±, whereas an on-shell particle is 
described by a positive p-spin only. 

The 3N basis states with one nucleon off mass shell 
were defined in Q. The case of the two-nucleon system 
with both nucleons off mass shell was discussed in [y]. 
Here we only need to combine the results of 0] and 0] in 
order to specify the more general case of 3N states with 
two nucleons off mass shell as required for the calculations 
of the 3N form factors. 

We use a basis of 3N states that are tensor products 
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of three one-nucleon states. The states are constructed 
in a sequence of steps. We start in the rest frame of the 
two-nucleon system composed of nucleons 2 and 3, and 
use a notation in which any variable with a "~" explicitly 
refers to this frame. Note that we use this notation only 
when the symbol without "~" is used to describe the same 
variable in a different frame. 

Nucleons 2 and 3 have four- momenta k 2 = (^20^2) 
and k 3 = (fc 30 ,k 3 ). Their relative momentum is 



P = \{k 2 ~ h) 



(23) 



and, since k2 + k 3 = 0, we have k2 = — k 3 = p. The 
energy components are subject to the constraint k 2 o + 
&30 = W, where W is the total energy in the two-nucleon 
rest frame. 

A particular state of nucleon 2, which will serve as our 
reference state, has its three-momentum p aligned along 
the positive z-axis, and its helicity and p-spin are A2 and 
p 2 . It is written as 



|(p,0,0)A 2 p 2 >, 



(24) 



where p is specified through its magnitude and polar an- 
gles 9 and 4> in the form (p, 9, </>). We use the same symbol 
for the four- vector and the magnitude of the three- vector 
in order to avoid an awkward notation. Which one is 
meant is always clear from the context. 

We denote a spatial rotation by an angle £ about an 
axis fi as i?(£n), and a boost along the direction h with 
rapidity 77 as B(r]h). The corresponding representations 
in Dirac space are 



S(R((n)) = exp -iS • h 



C 



. cos — — 1 sin , 

2 1 V 2 



!)E(ft-*) E ' 



(25) 



S(B(rjh)) = exp (^a ■ "■"J 

= lcoshm+sinhm£ 



n ■ e ) a; 



i=l 



(26) 



where er 1 , er 2 , er 3 are the Pauli matrices, and we introduce 
the Dirac matrices 









f 




0' 




Oil = 1 




0) 



The unit vectors 



(27) 



and e 3 point in the x, y, and z 



direction, respectively. The special combination of rota- 
tions R(ae 3 )R(f3e 2 )R(-fi 3 ) is abbreviated as Ra,0,-y 

With these conventions, a general state in which the 
three-momentum points in the direction of the polar an- 
gles 9 and cj) is 



|(p,M)A a pa) = S(R^ §0 )\(p,0,0)X 2 p 2 ) 



(28) 



The three-momentum of nucleon 3 points into the op- 
posite direction of the one of nucleon 2. In the phase 
convention by Wick [l2j, its state is defined by starting 
from the reference state (f2"4"]l with helicity A 3 and p-spin 
p 3 , applying a rotation Rn^,o, which aligns its momen- 
tum along the negative z-axis, followed by a rotation in 
a general direction, 



|(p,M)<W> 3 ) 



>S(R^ )\(p,n,ir)\ 3P3 ), (29) 



where s 3 = 1/2 is the spin of particle 3. The reason for 
introducing the extra phase factor e~ 17TSa is discussed in 
[H and Q 

The two-nucleon state can now be written as 

\(p,e,4>)X2X 3 ; P 2p S ) = e-^ S3 

xS(\j )0 ) {|(P, 0, 0), A2P2) \(P, 7T, tt), X 3P3 )} (30) 

where the common rotation g acts on the spaces of 
both particles simultaneously. 

Two-nucleon states with definite total angular momen- 
tum j and total helicity m arc obtained through 

\pjm,X 2 X 3 ;p2P3) = Vj I dd> I d9sm9 



xT>, 



0')* / 

m,A 2 — A 3 ' 



b,6,0)\(p,6,<j>)X 2 X 3 ;p 2 p 3 }, (31) 



where we use the familiar Wigner 2?-functions and the 
abbreviation 



Vj 



2j + 1 

47T 



1/2 



(32) 



We can treat these composite two-body states as if they 
belonged to a single particle with spin j and helicity m. 
Its total momentum is zero, because we are still in the 
two-body rest frame. Applying the boost Z(q) = B(i] q e 3 ) 
takes it to a state with three- momentum q in z-direction. 1 
The rapidity of the boost is determined by sinh?^ = 
q/W(q), where the invariant mass of the nucleon pair 



is W(q) = y[M t - E(q)] - q 2 . 

We proceed then by introducing the state of parti- 
cle 1, which in the three-body rest frame has its three- 
momentum in the opposite direction of the (23) pair, in 
exactly the same manner as the state of particle 3 was 
defined in the two-body rest frame, 



|(g,7r,7r)Ai;pi) 



1 5(i?^,o)|(g,0,0)A 1 ;p 1 ), (33) 



where the spin of particle 1 is si = 1/2. In this paper, 
nucleon 1 is always on mass shell. Its p-spin is therefore 
always pi = + and will be suppressed from here on. 



1 Note that this u q" is totally unrelated to the photon momentum 
"q" . It is unfortunate that both variables are traditionally given 
the same symbol, and we try to facilitate their distinction by 
using different font styles. 
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A tensor product of this state with the two-body state 
defined in Eq. (|3"Tj) yields a three-body state with zero 
total three-momentum and the momentum of the (23) 
pair, which has spin j and helicity m, in the positive 
z-direction, while nuclcon f has helicity Ai and its three- 
momentum in opposite direction to the pair. A general 
3N state is obtained by rotating the state aligned along 
the z-axis into a general direction by performing a rota- 
tion i?$,e.o- 

Repeating the procedure of Eq. (|3"Tj) , and exploiting 
the fact that a rotation about the z-axis commutes with 
a boost along the same axis, we can write the 3N partial 
wave helicity states with total angular momentum J and 
total helicity M in the form 

\qJM,pjm, X 1 (X 2 X 3 ); p 2 p 3 ) 

= rum J dS V%*_ Xl (S) £ d9 sin d { £ M _ X3 (0) 

xS(Rs) {\(q,ir,ir),\ 1 ) ® S(Z(q))S(R oAo ) 
x|(p,0,0),A 2 A 3;/92 p3)} (34) 
where we have used the abbreviations 



Rs — R$,e,< 



V 



(j). 



/•27T pTT /»27T 

dS = d<f> dO sin 6 / d<p . 
Jo Jo Jo 



(35) 



Note that so far we have suppressed the energy compo- 
nents g° and p a of the four- vectors, which don't play any 
role in the definition of the partial wave basis. While 
nucleon f is on mass shell and q° given by q° = E(q), 
p° is not determined through p when both pair nucle- 
ons are off mass shell, and p° has to be included in the 
specification of the state. 

Our isospin states are constructed as in 



[(t 2 t 3 )Tt!] TT Z ) , 



(36) 



where the isospins of nucleons 2 and 3 are first coupled 
to the pair isospin T, which is then in turn coupled with 
the isospin of nucleon 1 to yield the total 3N isospin T 
and its projection T z . Since for the 3N bound states 
t\ = t 2 = t 3 = T = 1/2, we will write the isospin states 
simply as \TT Z ). 

The complete specification of the helicity basis states 
is then 



\qJM,p pjm, Ai(A 2 A 3 ); p 2 pa;TT z ) 



rum 



dSV 



<M >m -xS S ) I ^sin^k-AaW 
Jo 

xS(R s ){\(q,Tr,K),\i)®S(Z(q))S(R ore!0 ) 

x\(p°p, 0, 0), A 2 A 3 ; p 2 p 3 )} <8> | [{hh)Tt x ] TT Z ) . (37) 

For the numerical solution of the Faddeev equation 
(jTTJ) for the 3N vertex function, we use a basis of states 



with good parity and particle exchange symmetry, which 
is obtained [3j through linear combinations of the states 
(Ell). 



III. RESULTS AND DISCUSSION 

In this section, we present the results of our calcula- 
tions of the elastic electromagnetic form factors of 3 He 
and 3 H in the complete impulse approximation (CIA). 
The underlying 3N vertex functions in (|22p were obtained 
by solving the corresponding Faddeev-type 3N bound 
state equation (fT"7|) which depends on NN scattering am- 
plitudes as dynamical input. The NN amplitudes are 
solutions of the CST two-body equation (jT8j) for given 
NN interaction models. We begin by a brief description 
of the interaction models used in our form factor calcu- 
lations. 



A. NN interaction models 

In order to study the model dependence of the form fac- 
tors, we used several members of a family of one-boson 
exchange (OBE) potentials that were first used in Ref. 
Q, in the first CST calculation of the triton binding 
energy. These models are based on the exchange of 6 
bosons, namely the tt, 77, a, 5, u>, and p. The free potential 
parameters were determined by fitting to the NN phase 
shifts below 350 MeV and to deuteron properties. 

A distinctive feature of these models is that they in- 
clude off-shell couplings of the scalar mesons a (isoscalar) 
and 8 (isovector) to nucleons, of the form 



g s A s (p',p) = g s 



1 



2m 



(38) 



where A s (p',p) is the meson- nucleon vertex, and p and 
p' are the four-momenta of the incoming and outgoing 
nucleons. The couplings proportional to v s do not con- 
tribute if the nucleons are on-shell, hence the name "off- 
shell couplings" . (Note that the definition of v s used in 
Eq. (|3"8"|) differs from that used in the recent work of Ref. 

It turned out that the triton binding energy is very 
sensitive to the strength of the scalar off-shell coupling. 
In the first exploratory calculations, it was convenient 
to vary v a and vs not independently, but keeping their 
ratio fixed. This was done by expressing each coupling 
constant through a common scale factor v, according to 



v a = —0.75 v , 



u s = 2.60 v. 



(39) 



The family of models discussed here has values of v vary- 
ing from to 2.6. 

Table U lists the parameters that determine the scalar 
off-shell couplings. The table also shows that, as v is 
increased from to 2.6, the triton binding energy varies 
over a large range. At the value v = 1.6 (model W16), 
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TABLE I: Scalar meson parameters of the NN potential mod- 
els used in the calculations of the 3N vertex functions. Also 
given are the xV^data of the models obtained in fits to a 
NN data base of 1994, and the corresponding triton binding 
energies E t . All masses and binding energies are in MeV. 





WOO 


W10 


W16 


W19 


W26 


V 


0.0 


1.0 


1.6 


1.9 


2.6 




5.84067 


5.50753 


4.99887 


4.67948 


4.05718 




0.0 


-0.75 


-1.2 


-1.425 


-1.95 




525 


515 


506 


501 


491 




0.14812 


0.69046 


0.62818 


0.47598 


0.25045 




0.0 


2.6 


4.16 


4.94 


6.76 


ms 


390 


540 


512 


474 


399 


xV N data 


3.00 


2.45 


2.25 


2.27 


2.56 


E t 


6.217 


7.411 


8.489 


9.072 


10.533 



the experimental binding energy is crossed, and simulta- 
neously the best fit to the NN data (lowest x 2 /Nd a ta) is 
achieved. Clearly, this singles out W16 as the most real- 
istic model. While xVNdata deteriorates when v moves 
away from 1.6, it does increase only moderately over 
the considerable range. We have therefore a convenient 
method at our disposal to generate a family of potentials 
that differ considerably in the 3N binding energy, but 
yield roughly equivalent fits to the NN data. This is an 
almost ideal situation to study the model dependence in 
our 3N electromagnetic form factor calculations. 

A more detailed description of these potentials and the 
complete list of parameters can be found in [l3| . In all of 
these models, the pion-nucleon coupling is of pure pseu- 
dovector form. (Recent fits to the data described in Ref. 
Q give x 2 /N — 1 and still show the same correlation be- 
tween Et and x 2 - Results from these new models will be 
discussed elsewhere.) 



B. Partial wave convergence 

The electromagnetic form factors of the 3N bound 
states have been calculated in the basis of partial wave 
helicity states described in Section III Dl With an in- 
creasing number of partial waves, the calculations be- 
come more and more time consuming. It is important to 
have an idea about the accuracy of the results that can 
be achieved with a limited number of partial waves. It 
is costumary to define a truncated basis of three-body 
states through the maximum value of the two-body total 
angular momentum, j max , of the included basis states. 
Table |n] shows the charge form factor of 3 He, for a par- 
ticular model calculation which will be described later, at 
5 different values of the momentum transfer Q, depend- 
ing on j max . Table ITTT1 shows the same for the magnetic 
form factor of 3 He. 

What both tables indicate is that at lower Q a good 
stability of the results is reached already with j max = 2. 



TABLE II: Partial wave convergence of the charge form factor 
of 3 He at selected values of the momentum transfer Q. The 
first column shows the maximum total pair angular momen- 
tum included in the 3N partial waves. Columns two to six 
are the form factor calculated in the corresponding truncated 
partial wave basis at Q=1.0, 3.0, 5.0, 7.0, and 9.0 fin" 1 . Ex- 
cept for those in column two, the values of the form factor 
are multiplied by a power of ten indicated in the last line. 
All results are for model W16 with the MMD nucleon form 
factor. 



Jmax 






Q (mr" 1 ) 








1.0 


3.0 


5.0 


7.0 


9.0 



1 


0.5585 


1.7603 


-1.5546 


-0.9903 


1.039 


2 


0.5716 


1.8772 


-1.5314 


-0.9815 


1.201 


3 


0.5733 


1.8866 


-1.5519 


-1.0471 


1.174 


4 


0.5744 


1.8903 


-1.5511 


-1.0618 


1.119 


5 


0.5743 


1.8904 


-1.5515 


-1.0623 


1.128 


6 


0.5744 


1.8896 


-1.5518 


-1.0644 


1.119 






xl0~ 2 


xlO" 3 


xl0~ 4 


xlO" 5 



TABLE III: Same as TableHU but for the magnetic form factor 
of 3 He. 



'max 






Q (mr" 1 ) 








1.0 


3.0 


5.0 


7.0 


9.0 


l 


0.5063 


-0.6586 


-3.4357 


-0.6531 


4.3456 


2 


0.5202 


-0.6207 


-3.6999 


-1.0678 


4.4701 


3 


0.5217 


-0.6181 


-3.6482 


-0.8808 


4.7954 


4 


0.5227 


-0.6200 


-3.6497 


-0.9454 


4.5796 


5 


0.5226 


-0.6190 


-3.6489 


-0.9376 


4.6228 


6 


0.5226 


-0.6206 


-3.6479 


-0.9413 


4.5910 






xlO" 2 


xlO" 3 


xlO" 4 


xlO" 5 



Not surprisingly, the results are less stable for higher Q, 
although uncertainties of the order of a few percent are 
perfectly satisfactory. 

Diagrams A and the sum of B and C require numerical 
integrations over two momenta and two angles, while in 
diagrams D and the sum of E and F two momenta and 
three angles are integrated. In principle, the precision of 
the numerical integrations can be improved upon by in- 
creasing the number of integration points. However, we 
feel that the current precision, which we estimate to be 
about three significant figures, is sufficient for the pur- 
pose of this work. Small changes of the order of one 
percent would not be visible in the figures presented. 

C. Model dependence of the CIA 

Figures SH7] compare our CIA results for the selected 
models listed in Table U to IARC (Impulse Approxima- 
tion with Relativistic Corrections) calculations within 
the framework described in Refs. [1, [Hj]. The IARC cal- 
culations use a one-nucleon-current with wave functions 
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FIG. 4: (Color online) Charge form factors of 3 He and 3 H for different NN interaction models in CIA. The left panel shows 
the form factors divided by a scale function (solid straight line in the right panel). Shown are the results for models WOO 
(short dash-dotted), W10 (long dash-dotted), W16 (dash-double dotted), W19 (long dashed), and W26 (short dashed), all with 
the MMD nucleon form factor. In addition, W16 (solid line) is compared to IARC by Marcucci (double dash-dotted), where 
both calculations used the Galster nucleon form factor. The theoretical error bars of the Greens function Monte Carlo IARC 
calculations are also given. The right panel shows the W16/Galster and IARC/Galster results together with the scale function 
(solid straight line) in the traditional semi-log plot. Both panels also show the experimental data (full circles) [3] . 



obtained from the Argonne AV18 NN and Urbana IX 3N 
potentials, and include first-order relativistic corrections. 
Our models do not include the strong Coulomb correc- 
tions for the pp interaction, and hence give the same bind- 
ing energies for both 3 H and 3 He. To help in the compar- 
ison of our results with the IARC, L. Marcucci provided 
us with IARC calculations in which the pp Coulomb in- 
teraction was excluded [l6| . such that the used 3 H and 
3 He wave functions also had the same binding energies. 

Figures 2] and [3] show the charge and magnetic form 
factors of 3 H and 3 He, while Figs. [5] and [7] show the 
isoscalar and isovector combinations of the charge and 
magnetic form factors. Each figure has four panels. The 
right hand panels give a log plot of the absolute value 
of the form factors vs. Q. Note that they all fall off by 
about 6 orders of magnitude over the range in Q shown 
in the plot. In each of these right hand panels an ex- 
ponential scale function is shown (the straight line). To 
get a better idea of the relative differences in the models, 



the left hand panels show the form factors divided by the 
scale function. This removes the strong exponential de- 
pendence and permits us to show both the sign and the 
relative value of each curve on a linear scale. 

The first conclusion that can be drawn from these 
graphs is that, in every case, model W16 is remarkably 
close to the IARC calculation at low Q, confirming that 
these two calculations are in essential agreement with one 
another. This is a very pleasing result: it confirms for 
the first time that the CST not only yields the correct 
3N binding energy, but that it is also capable of a good 
description of the electromagnetic structure of the 3N 
bound states. The differences in the underlying dynam- 
ics between IARC (the nonrelativistic AV18 NN potential 
and an irreducible 3N force) and the CST (relativistic NN 
kernel with off-shell couplings, no irreducible 3N force) 
seem to be less important at small Q than the fact that 
they yield the same nuclear binding energy. As will be 
discussed in more detail in Sec. IIII F[ in the case of the 
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FIG. 5: (Color online) Magnetic form factors of 3 He and 3 H for different NN interaction models in CIA. The various curves 
are defined in the caption of Fig. [4] 



family of CST interaction models used here, the CST 
single-nucleon current contains approximately the same 
physics as the one of I ARC. The close agreement of the 
results is therefore perfectly understandable. 

Next, note that model W16 departs from the IARC at 
a Q of about 4 to 6 fm _1 . The Monte Carlo calculations 
of the IARC begin to show large errors for Q > 6 fm _1 , 
so it is difficult to make a precise comparison above this 
point, but the trends are still clear. The charge form fac- 
tors agree to Q ~ 6 fin -1 beyond which the IARC mod- 
els oscillate with a shorter wavelength than does W16. 
A similar thing happens with the magnetic form factors, 
but the IARC and W16 break away at a lower Q ~ 4 
fm _1 . The isoscalar and isovector combinations show 
the same behavior. 

The variation in the CIA models is remarkably smooth, 
and seems to be explained almost entirely by the different 
three-body binding energies predicted by each of these 
models. Recall that the low momentum behavior of the 
charge form factors can be written 



Fc(Q 2 



1 



■Q 2 (r 2 ) + 



(40) 



Using the values of the root mean square charge radii 
in Table IIV1 which also shows the magnetic radii for 



completeness, it can be verified directly in Fig. [5] that 
the well-known relation (V 2 ) cx 1/E t between the mean 
square radius of the bound state and the binding energy 
E t holds very well for the CST models. 

This behavior predicts model WOO with the smallest 
binding energy will have the greatest curvature at Q = 0, 
and that W26 with the largest binding energy will have 
the smallest curvature, a behavior confirmed by the plots. 
The surprising fact is that this behavior seems to persist 
on to higher Q, with WOO oscillating with the shortest 
wave length and W26 the longest. 

Finally, note that the comparison with the data shows 
significant discrepancies beyond Q ~ 1 fm _1 , particu- 
larly for the isovector magnetic combination. In the case 
of the IARC calculations, it has been shown that the 
inclusion of large exchange-current contributions brings 
the theoretical calculations into good agreement with the 
data. While we are not yet ready to calculate the inter- 
action currents for our CST models, we believe that they 
will be large also. We reemphasize that the observed dis- 
crepancy between the data and the theoretical results is 
due to known physics not included in the impulse approx- 
imations in both cases. This aspect is discussed further 
below in Sec. ME 
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FIG. 6: (Color online) Isoscalar and isovector charge form factors of the 3N bound states for different NN interaction models 
in CIA. The various curves are defined in the caption of Fig. [3] 



TABLE IV: Charge and magnetic root mean square radii of 
3 He and 3 H for the five models of Table [TJ calculated in CIA 
with the MMD parametrization of the nucleon electromag- 
netic form factors. The corresponding 3N binding energies 
are also listed. 



r.m.s. radius (fm) 
charge magnetic 



Model E 


t (MeV) 


3 He 


3 H 


3 He 


3 H 


WOO 


6.217 


2.204 


1.986 


2.361 


2.172 


W10 


7.411 


2.030 


1.842 


2.195 


2.027 


W16 


8.489 


1.917 


1.742 


2.061 


1.946 


W19 


9.072 


1.869 


1.698 


1.991 


1.914 


W26 


10.533 


1.783 


1.613 


1.819 


1.862 



TABLE V: Magnetic moments of 3 He, 3 H, as well as their 
isoscalar and isovector combinations and fiv in nuclear 
magnetons. The first two lines are calculated with model 
W16 in CIA, in combination with the MMD and the Galster 
parametrization of the nucleon electromagnetic form factors. 
The third line is a IARC/Galster calculation by Marcucci. 
The last line shows the experimental values. 





M 3 He) 


m( 3 h) 


MS 


fJ-v 


W16/MMD 


-1.747 


2.550 


0.402 


-2.149 


W 16/ Galster 


-1.749 


2.546 


0.398 


-2.147 


IARC/Galster 


-1.763 


2.572 


0.404 


-2.168 


Experiment 


-2.127 


2.979 


0.426 


-2.553 



D. Contributions of the six diagrams 



Finally, in Table|V]we compare calculations of the mag- 
netic moments of the three-body nuclei. Note the reason- 
ably close agreement between the calculations, and the 
well known discrepancy with the data (due to missing 
exchange current contributions). 



The contributions of each of the six diagrams that 
make up the CIA are shown in Fig. [5J Note that all 
give comparable results; none can be neglected. 

It is particularly interesting to see that the sum of di- 
agrams B and C (each contains an integrable singularity, 
but their sum has no singularity) is almost precisely equal 
to diagram A (in all cases) . This was expected, but a sim- 
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FIG. 8: (Color online) Mean square charge radii of 3 He and 
3 H versus inverse binding energy for models WOO, W10, W16, 
W19, and W26, in order from right to left. The straight lines 
are linear fits through the points and show that the expected 
linear relation between (r 2 ) and 1/Et is indeed well satisfied 
for these models. 



ilar result does not hold for the sum of diagrams E and 
F. Each of these is singular, their sum is finite, but it is 
not equal to diagram D, as originally expected. Note also 
that the contributions of diagrams A and D to the charge 



form factors are very similar, at least at low Q 2 , but that 
they give substantially different results for the magnetic 
form factors. These relationships will be studied in a 
future paper. 



E. Dependence on the nucleon form factors 

In Ref . it was shown how to construct gauge invari- 
ant interactions using composite particles. Subsequently 
[l7l HH the general form of the one nucleon current op- 
erator for use in such a program was derived. The one 
used most frequently is 

f N (p',p) = fo(p' 2 , P 2 ){(Fi(Q 2 ) i)t + r} 

+go(p' 2 ,p 2 )A-(p'){(F 3 (Q 2 )~l)r + r}A-(p) 

+f^ 2 , P 2 )F 2 (Q 2 )^p^, (41) 
2m 

where Fi^(Q 2 ) are the on-shell nucleon form factors, 
F^iQ 2 ) is a completely unknown form factor describing 
the off-shell structure of the nucleon (subject to the con- 
straint that F 3 (0) = 1), A_(p) = (m-f)/2m, 7^ = 
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FIG. 9: (Color online) Contributions of the various diagrams to the charge (upper panel) and magnetic (lower panel) form 
factors of 3 He and 3 H. The solid line is the total result in CIA (sum of diagrams A - F) for the NN interaction model W16 
and the MMD nucleon electromagnetic form factors. The other lines are the partial results for diagrams A (dashed), B+C 
(dash-dotted), D (double dot-dashed), and E+F (double dash-dotted). 



7 M — q^tji/q 2 , /o is a completely undetermined function 
of p' 2 and p 2 (subject to the constraint /g(m 2 , to 2 ) = 1), 
and fo and go are functions of p 2 and p' 2 completely de- 
termined by the Ward-Takahashi identity for a dressed 
nucleon 

WW,P) = S^(p) - fi^V) (42) 
where Sn(p) is the dressed propagator for the nucleon 

Hp 2 )] 2 



S N {p) = 



(43) 



In the applications discussed in this paper, the nucleon 
form factor is 



Kp 2 ) = 



(A 



m 



2\2 



(A? 



l n — to 2 ) 2 + (m 2 — p 2 ) 2 
and /o and go are given in terms of h 



(44) 



h to 2 



■p 



9 (P ,P)= 77 T 



12 h! TO 2 



V 



h p' 2 — p 2 
4m 2 



p' 2 — p 2 



(45) 



with h = h(p 2 ) and hi = h(p' 2 ). 

In order to include the isospin dependence of the cur- 
rent correctly, care has to be taken to maintain the Ward- 
Takahashi identity satisfied. The following expression 
meets this requirement: 



3n(p'>p) = 



1 



(F lp (Q 2 

F 2p (Q 
1 + r 3 



-fo(p' 2 , P 2 h 
1 + r 



+ F ln (Q 2 ) 



1-T 3 



h(p'\p 2 )T 



1 + T A 



F 2n {Q 2 )- 



g (p' 2 ,p 2 )A_(p>h»A_(p) 



(F 3p (Q 2 ) - 1) 



T " - ' o2\ ^ 



+ F 3n (Q 2 )- 



xg (p' 2 ,p 2 )A-(p r )rA-(p) 



(46) 



Note that the nucleon current is always used in con- 
junction with a conserved electron current, and hence 
the terms in g M vanish. Making this simplification in 
Eq. (|46p . and combining the proton and neutron current 
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FIG. 10: (Color online) Charge (upper panels) and magnetic (lower panels) form factors of 3 He and 3 H in CIA, divided 
by the respective scale functions of Fig. 2] with model W16 and different nucleon electromagnetic form factors. The MMD 
parametrization is used for the on-shell part of the nucleon current in all cases. The solid line is obtained using only the on-shell 
nucleon current NCI (i.e., fo = fb = 1 and go = 0). The dashed line uses the full current (|4ip in the version NCII (with fb = fo 
and i*3jv = Gen)- Finally, the dash-dotted line shows the full current NCIV (with fb = 1). We also tried replacing F$n — Gen 
by F3N = -Fijv for each case (currents NCIII and NCV); the effect is too small to be seen on the plot. 



into a single, isospin dependent expression, gives 

+ r o{p '\^)F 2N {Q 2 ) 1 -^ 

+go(p' 2 ,p 2 )F 3N (Q 2 )A^p>)rA„(p) , (47) 

where, for i = {1, 2, 3} and the nucleon isospin projection 
r 3 , 

F lN {Q 2 ) = F ip {Q 2 ) 1 -^- + Fm (Q 2 )^—^ ■ (48) 

In previous applications, fb has been taken to be equal 
to fo, and usually F 3 n — Gen (where N stands for p or 
n), but no systematic study of the dependence of these 
factors has been completed. 

In this work, we performed 3N form factor calcula- 
tions with 5 different single-nucleon currents, which can 
be characterized by the different choices of the off-shell 
form factors fo, fb, go, a s well as of i^jv- 



The first nucleon current, labeled NCI, is the usual on- 
shell current, with /o = /o = 1 an< i So — 0. Next we con- 
sider a genuine full off-shell current, NCII, with fo and 
go given by Eq. gS]), fb = fo, and F 3N = G E n- Current 
NCIII differs from NCII only by the choice F 3N = F 1N . 
To study the sensitivity of the results on fb we con- 
structed current NCIV, which is equal to NCII except 
for fb = 1. Finally, current NCV equals NCIV apart 
from F 3 n = Fin. 

Figure [TU] shows the effect of some variations in the 
choice of the single-nucleon form factors. In general, the 
changes are very small for Q < 5 fm . Only at larger 
values of Q the various cases begin to diverge from each 
other, particularly in the charge form factors. The three 
displayed lines correspond to NCI, NCII, and NCIV. The 
results for NCIII are indistinguishable from NCII in the 
figure, as are those for NCV from NCIV. We conclude 
that the inclusion of off-shell form factors, as well as the 
choice of fb(p' 2 ,P 2 ), has only marginal effects at low Q, 
but becomes relevant for larger values of Q. On the other 
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FIG. 11: (Color online) Examples of Z-graphs in time- 
ordered perturbation theory, which are automatically in- 
cluded through the coupling of the photon to negative-energy 
nucleon states in the full CIA calculations. In the 3N system, 
the meson lines (dashed) can couple to any of the remaining 
two nucleons (not shown). 




(b) 

FIG. 12: (Color online) The pion interaction current. With 
ps coupling the contact current (b) is zero and the Z-diagram 
(a) is large. With pv coupling, the Z-diagram (a) is small, 
and the contact term (b) is large. There is an approximate 
equivalence between (a) calculated using ps coupling and (b) 
calculated with pv coupling. 



hand, the particular parametrization of i^3jv(Q) is much 
less important. 



F. Negative energy states as relativistic corrections 
vs. interaction currents 

We have seen that the CIA models and the nonrela- 
tivistic I ARC agree quite well at low Q, so that (i) dif- 
ferences in the physics they describe shows up only at 
higher Q, and (ii) their failure to explain the low Q data 
is probably due to the omission of the same physics. Here 
we will look at both of these questions. 

One must be cautious in comparing these two ap- 
proaches - the language used to describe the physics dif- 
fers and the same processes can be called by two dif- 
ferent names. Negative energy states are automatically 
included in the CST calculations, and are considered rel- 
ativistic corrections, but in the IARC formalism they are 
considered to be interaction currents. 

The negative energy (or interaction current) terms we 
are discussing are illustrated in Fig.QTj The figure shows 
a selected set of Z-diagrams involving one, two and three 
pion exchanges. These are "time-ordered" diagrams with 
time flowing from left to right, so the diagonal lines (flow- 
ing backward in time) represent negative energy contri- 
butions, which can be reinterpreted as contributions from 
the production of virtual nucleon-anticuclean pairs. In 
the Feynman diagram formalism, more natural to the 
CST, these Z-diagrams are automatically included as 
part of the off-shell nucleon propagator. The CIA au- 
tomatically includes contributions from the infinite sum 
of all of these diagrams. In the IARC formalism, only 
diagram [TIT a) is included. The other processes, thought 
to be small, are never included in IARC. Furthermore, 
in the nonrelativistic language of IARC, some of these 
higher order diagrams generate three-body currents, so, 
in the language of IARC, our calculation also includes 
some three-body currents. However, in the language of 
the CST, none of these contributions are three-body cur- 
rents. 

There is another ambiguity that adds to the confusion. 
If pseudoscalar (ps) pion coupling is used to calculate the 



pion exchange currents (the choice made by IARC) , then 
there is no four-point jttNN current, but when pseu- 
dovector (pv) coupling is used, there is a (large) four- 
point jttNN current. This is illustrated in Fig. [12J It 
has been known for a long time that the leading order 
contribution from the Z-diagram (a) in ps coupling is 
equivalent to the contact current (b) in pv coupling. In 
the CST models used here, a pure pv pion coupling was 
used, so the Z-diagrams are small and the interaction 
current of Fig. \V2\b) must be added (it would be part of 
the contributions coming from Figs.[3]GI - J). Hence, our 
results are similar to the IARC calculation because, even 
though we include the Z-diagram |T2](a) and IARC do 
not, ours is small. The large part of the pion interaction 
current, [12(a) for IARC and HHb) for CST, is omitted 
by both calculations. 

How big are the negative energy contributions shown 
in Fig. [TTI ? We have said that they are small, because 
we use pv coupling for the pion. A more precise picture 
is given Fig. Q21 which shows the influence of negative- 
energy states in three-body channels on the charge and 
magnetic form factors of 3 He and 3 H. It compares the full 
calculation with W16/MMD with another calculation in 
which only positive-energy three-body channels are taken 
into account. Note that this does not mean that negative- 
energy states are completely eliminated from these calcu- 
lations, because the three-nucleon vertex functions were 
calculated from all channels. But this comparison does 
show the effect of removing all of the diagrams in Fig. 1111 
The figure shows that the charge form factors of 3 He and 
3 H are slightly changed, mainly for momentum transfer 
larger than about 6 fm _1 , while the magnetic form fac- 
tors remain essentially unchanged. 



G. Critique of the physics missing from the CIA 

We conclude with a brief discussion of the physics miss- 
ing from the CIA. It is useful to look at the isoscalar and 
isovector magnetic and charge combinations separately. 

Careful examination of the isoscalar magnetic form 
factor given in Fig. [7] shows that the CIA is very close to 
the data, except possibly in the region near Q from about 
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FIG. 13: (Color online) Charge (upper panels) and magnetic (lower panels) form factors of 3 He and 3 H divided by the respective 
scale functions of Fig. [3] The solid line is the full result, the dashed line excludes all three-body negative-energy states in the 
3N vertex function. The results are obtained with the NN model W16 and the MMD nucleon form factor in CIA. 



2 to 3 fm -1 (where the discrepancy is about ~ 20%). 
This is confirmed by the IARC studies 0, which show 
that interaction currents in this channel are quite small 
(the largest "corrections" to the one-body current are the 
spin-orbit relativistic corrections, included in both the 
IARC and in the CST models). The model dependence 
shown in the figure suggests that the discrepancy in the 
region Q from 2 to 3 fm -1 (not visible on the logarithmic 
scale shown in the right hand panel and in Ref. |9|) might 
very well be corrected by use of more accurate wave func- 
tions. Accurate 3 H wave functions were produced after 
this work was finished, and a better 3 He wave function 
(with the correct binding energy) can be produced once 
we develop an accurate pp interaction model. 

The isovector magnetic form factor (Fig. [7]) and both 
the isoscalar and isovector charge form factors (Fig. O 
show large discrepancies that are explained in Ref. [9| by 
large pion and rho interaction currents missing from both 
the IARC and the CIA results presented here. It remains 
to be seen whether these corrections, when added to the 
CIA, will bring the CST predictions into agreement with 
the data. 

Finally, as observed in all of the figures above, rela- 



tivistic effects grow significantly with increasing Q 2 , and 
the new Jefferson Laboratory high Q 2 measurements of 
the 3 He form factors, to be released soon, will provide 
a strong test of the relativistic theory. The existence or 
non-existence of a second minimum will be a particularly 
interesting signal. Our CIA calculations show a second 
minimum, but until we have added the missing interac- 
tion currents we are unable to predict its location. 



IV. CONCLUSIONS 

We have performed the first numerical calculations of 
the electromagnetic three-nucleon form factors in the 
Covariant Spectator Theory. This framework is man- 
ifestly covariant and includes relativistic effects such 
as boosts, Wigner spin rotations, and negative-energy 
states exactly, without resorting to expansions in orders 
of v/c. These calculations were done in Complete Im- 
pulse Approximation for a number of relativistic three- 
nucleon bound-state vertex functions, derived from a 
family of two-nucleon interaction models that produce 
different three-nucleon binding energies while maintain- 
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ing a roughly equivalent good fit to the nucleon-nucleon 
data. From the findings of these calculations, we want to 
highlight the following: 

The first observation is that the CST description of the 
three-nucleon electromagnetic form factors yields very 
reasonable results, confirming the validity of this formal- 
ism. The calculations are demanding, but numerically 
stable. 

Our study of the model dependence of the results indi- 
cates that the three-nucleon binding energy determines 
their behavior to larger momentum transfer than previ- 
ously expected. This is also evident when we compare to 
nonrelativistic impulse approximation calculations with 
relativistic corrections (IARC), which are based on the 
Argonne AV18/UIX two- and three-nucleon potentials. 
The three-nucleon form factors obtained with model 
W16, which reproduces the experimental triton binding 
energy and comes also closest to the one predicted by 
the nonrelativistic AV18/UIX potentials, are remarkably 
similar to the IARC form factors up to moderate values 
of the momentum transfer (around 4 to 6 fm -1 ). 

We find that each of the six diagrams of CIA is signif- 
icant, so none can be omitted. 

Variations in the parametrization of the single-nucleon 
electromagnetic form factors, as well as the inclusion of 
off-shell form factors, have little effect on the results. 

A similar observation can be made about the influence 
of negative-energy states in the three-nucleon channels, 
that give rise to Z-graph-type contributions to the form 
factors. We explain their smallness by the use of pure 
pseudo-vector ttNN coupling in our two-nucleon poten- 
tials, which suppresses Z-graphs. The same reason makes 
it perfectly understandable why our results are so similar 
to the ones in IARC, that do not include any Z-graphs. 

Since the CIA calculations do not contain interaction 
currents, which we believe to be large, a good agreement 
with the data over the whole range of the considered mo- 
mentum transfer cannot yet be expected. An execption 
is the isoscalar magnetic form factor, where interaction 
currents are small, and indeed we find a good agreement 
of our results with the data. In all other cases, interac- 
tion currents will have to be calculated when a very close 
description of the data for all values of Q is the objective. 
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APPENDIX A: CALCULATION OF THE 3N 
VERTEX FUNCTION WITH TWO OFF-SHELL 
NUCLEONS 

As discussed above, the contributions when two nucle- 
ons are off-shell require knowledge of the two-body off- 
shell vertex function defined in (|2H)) . This is determined 
by quadratures from the two-particle off-shell scattering 
amplitude, which is in turn determined by quadratures 
from the two-body off-shell kernel Vp a ^ a i(k 2 , k' 2 ;P23), 
as illustrated in Fig. O 

The off-shell kernel enlarges the number of degrees of 
freedom in two different ways. 

First, since the kernel connects to the off-shell propa- 
gator for particle 2, its expansion [similar to the expan- 
sion (|15p for particle 3] requires knowledge of both the 
positive and negative energy projections of the kernel: 



u P 2 (k 2 ,\ 2 )V^y 2a ,(k 2 ,k' 2 ;P 23 ). (Al) 



The p 2 = — projections, not needed before, are defined 
in Ref . Q , and can be computed straightforwardly from 
knowledge of the Dirac structure of the OBE model. 

Second, the relative energy of the final state is now no 
longer fixed. In general, this relative energy is 



^23 



2 (*2 



*8) 



\W 2 . 



(A2) 



When particle 2 is on shell, k 2 = Ek 2 and the relative en- 
ergy depends on | ki2 1 and is not an independent variable. 
When both particles are off-shell, k 23 is an independent 
variable. In this work we found it convenient to write k 23 
in terms of the new variable Xq, where 



k° 

^23 



Xq 



[Eh 



\W 2 ; 



(A3) 



so that xq = 1 when particle 2 is on shell. 

To have the behavior required by the (generalized) 
Pauli principal, the potential must exhibit a Pauli ex- 
change symmetry (meaning that, depending on the 
isospin, it must be either symmetric or antisymmetric) 
when particles 2 and 3 are exchanged. However, when 
particles 2 and 3 are exchanged, k 23 — > —k 23 , so that 
xq — > — xq- We conclude that the kernel must have the 
Pauli exchange symmetry when xq — > — xq. For the on- 
shell case (xg = 1), this symmetry is discussed in detail 
in Refs. HQ. 

This symmetry is imposed in practice by explicitly 
symmetrizing the OBE kernels. The symmetrized me- 
son propagators are 



1 



± 



2 \m% - t+(xo) ml ~ t-(xo) 



(A4) 



where q± = ±k 23 — k' 23 and t± = q±- In Refs. 0, Q the 
first term in (| A4[l is referred to as the direct term; the 
second as the exchange or "alternating" term. Since the 
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initial particle 2 is on-shell, the momentum transfers are 



t±(x ) 



E' -\W 23 TXo{E-\Wm) 
-(kTk') 2 , 



(A5) 



where, for simplicity, we use the notation k = k 2 % and 



k' = 
W 23 : 



v 23- 



Near the physical scattering region, when 



2E ~ 2E', the momentum transfers t±(xo) 



— (k =F k') are negative. However, when the momenta 
are very large, t± can become positive, and the denomi- 
nators rr? v — t± can be zero. To see that the momentum 
transfers can become positive, consider the case when 
k - k' = and Ikl -> oo. Then 



lim t±{x ) -> k 2 (l ±x ) 5 



(A6) 



the on-shell amplitudes. This prescription can be sum- 
marized by the replacements t±(xo) — > i±(±l) which 
insures that the meson propagators for both the direct 
(same particle on-shell) and "alternating" (different par- 
ticles on-shell) have no singularities. The exact xq depen- 
dence is retained in the numerators N(q) so that some 
differences between direct and alternating terms are pre- 
served. Hence, under prescription A all symmetries are 
preserved, but some terms (those that would arise from 
antisymmctrizing the denominators) are zero. In choos- 
ing prescription A we lose some of the off-shell depen- 
dence. A more complete picture can be obtained from 
prescription C (described in Ref. @), and the conse- 
quences of this choice will be described elsewhere. 



Hence the momentum transfer always becomes positive, 
and the kernel has a singularity, unless xq — 1 and the 
same particle is on-shell both before and after the inter- 
action (the direct term with t = t + ). In this case it is 
easy to show the t+(l) < for all k and k'. 

A detailed analysis shows that some of these singu- 
larities arise from physical particle production; others 
are spurious singularities that are cancelled by higher 
order kernels neglected in the OBE approximation. In 
either case these singularities should be removed, and 
prescriptions for removing them are discussed in detail 
elsewhere (the interested reader should study Ref. [(|). 
The work presented in this paper uses prescription A for 



APPENDIX B: CALCULATION OF THE 
THREE-NUCLEON ELECTROMAGNETIC 
CURRENT IN CIA 

In this appendix we show some details of the calculation 
of the 3N electromagnetic current in CIA, i.e., of dia- 
grams (A) to (F) of Fig. [3l The calculations are carried 
out in the Lab frame. Of course, since we are working in 
a covariant framework, the final result does not depend 
on this particular choice of reference frame. 

First, we write the helicity partial wave basis states 
(|3T[) in a way that shows the involved transformations on 
the individual particle states in more detail: 



J 



qpPpM jm\i\ 2 \ 3 p 2 p 3 TT z ) 



2j + l 



d<P / dO sin O 
lo Jo 



d6 sin ev 



(1/2)* 
M,m— A 



\pi) <8> S(R*,e,o)S{RK,n,o)u + {q, Ax) <8> |p 2 > ® S{R w )S(Z{q))S{R i§fi )u p2 (p, A 2 ) 



\ps) ® S(B*,e, )S(Z(q))S(Ri s )S(i?^, K 3 (p, A 3 ) 



\T%) 



(Bl) 



The total angular momentum is always J = 1/2 for 
the 3N bound states and is from here on suppressed 
in the state kets. The possible values of the other dis- 
crete quantum numbers are M = ±1/2, j = 0,1,..., 
m = + Ai, A 2 , A 3 = ±1/2, subject to 

the constraints \m — Ai| < 1/2, |A 2 — A3] < j, and p 2 , 
P3 = ±. The four-momenta pi, p 2 , p 3 of the three nucle- 
ons are obtained by the space-time transformations 

Pi = R*,e,oR*,*,o(E(q),Q,0,q) , (B2) 

P2 = R^, ,oZ(q)R < j > g Q (p°,O,O,p) , (B3) 

p 3 = R^0, o Z(q)R lSo R w ^ o (W(q) -p°,0,0,p) , (B4) 

and for single-particle state vectors we use the normal- 
ization (p'|p) = (2tt) 4 ^ 4 (p' -p). 



The helicity spinors in (|B1[) are 

u + (p,X) = I n "TJ- I ®x(A). 



cosh ^ 
2A sinh \ 



u (p, A) = 



-2Asinh^ \ ... 
cosh f 2 )^ (A) ' 



(B5) 



with the rapidity rj p given by tanhry p = p/E(p), and the 
two-component spinors 



X(V2) = 



x("l/2) = 



(B6) 



The state (|B1[) has total angular momentum 1/2 and 
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helicity M. As long as the 3N bound state remains at 
rest or is boosted along the positive z-axis, M is also the 
z-projection of the total angular momentum. 

The total four-momentum of the 3N bound state at 



rest is Pt = (M t , 0,0,0). Its vertex state with nucleon 
1 on mass shell, helicity M, and isospin projection T z , 
can be written as a linear combination of the basis states 
(El): 



x CiqfpMjmXxMXz&p'zTT^Op^ (p, A a )O wp j (p, X 3 )\qp°pMjmX l X 2 X 3 p 2 p 3 TT z ) , (B7) 

I 




where the functions C are real because of time reversal 
invariance, and the matrix O is defined as 



CV P (p,A)=^ (p, A)«"(p, A) 

\ ~Z/\- —1 / 

\ m / p'p 

(B8) 

The integration over </ in (1B7|) extends up to the finite 
value q s = [M 2 — m 2 )/2M t for which the mass of the 
2N subsystem becomes zero and the limit of the region 
of timelike 2N states is reached. As was shown in 
the functions C go to zero smoothly as q approaches q s , 
which makes it possible to treat q s as a natural cut-off 
momentum without making the 3N vertex functions dis- 
continuous. 

When the Feynman diagrams of the CIA are evalu- 
ated, one encounters frequently the particular sequence 
of Lorentz transformations of a boost in z-direction, fol- 
lowed by a rotation about the y-axis and then by another 
boost in z-direction. It turns out to be very useful that 
this is equivalent to one boost in z-direction between two 
rotations about the y-axis, in the following way: 

Bfae^JZC^Bfoe 3 ) = R(6 2 e 2 )B( V3 e 3 )R(6 3 e 2 ) . 

(B9) 

For given rapidities 771, 772, and rotation angle 6±, the 
corresponding rapidity 773 and rotation angles 8 2 and #3 
can be found from: 

sinh 773 = \/ (cosh 771 cosh 772 + sinh 771 sinh 772 cos #1 ) 2 — 1 . 

(BIO) 

If 773 > then 



with < 6» 2 , 9 3 < 2tt. If 773 = then 6 3 = and 



sin t)2 



COS O2 



sin t> 3 



COS #3 = 



sin 0\ sinh 772 
sinh 773 

sinh 771 cosh 772 + cosh 771 sinh 772 cos 6± 
sinh 773 

sin 0i sinh 771 

sinh 773 ' 
cosh 771 sinh 772 + sinh 771 cosh 772 cos 0\ 
sinh 773 



(BH) 
(B12) 
(B13) 
(B14) 



sin O2 = sin V\ , 

COS 02 — COS 0i 



(B15) 
(B16) 



where < 2 < 2n. 

Similarly useful, a rotation about the z-axis between 
two rotations about the 7/-axis can be replaced by a ro- 
tation about the y-axis between two rotations about the 
z-axis: 

R(8 1 e 2 )R((f> 1 e 3 )R{02e 2 ) = R(0 2 e 3 )R{0 3 e 2 )R{cj) 3 e 3 ) . 

(B17) 

From the given angles 0\ 1 02, and 4>i one obtains first 

6>3 = arccos(cos 0\ cos #2 — sin#i sin 02 cos^i) . (B18) 
If 3 ^ and 3 then 
sin 8 2 sin <f>i 



sin <p2 

COS (f>2 

sin cj> 3 
cos (j) 3 



sin v 3 

sin 8\ cos 02 + cos 0\ sin 02 cos 6\ 



sin 6*3 



sin c/i sin 0i 



sin v 3 

cos 0i sin 02 + sin 0\ cos 02 cos cj)\ 
sin 03 



(B19) 
(B20) 
(B21) 
(B22) 



with < <j>2, 4>3 < 27r. On the other hand, if cos #3 = ±1 
then 03 = 0, and 



sin <p2 = - sin (f>i (cos 8\ ± cos 2 ) ■ 



(B23) 



cos 4>2 = — [cos 0i (1 ± cos 0i cos 8 2 ) T sin 6\ sin 8 2 ] , 

(B24) 

with < 02 < 2n. 

The 3N vertex functions were calculated in a helicity 
partial wave basis, which means that we have to calculate 
the matrix elements of the electromagnetic 3N current in 
terms of the functions 

C(qp°pM jm\i\2\ 3 p2P3TT z ) = 

(qp°pMjm\ 1 \ 2 \ 3P2 p 3 TT z \r 1 (P t ,M,T z )) . (B25) 
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Formally, we can write each of the diagrams of the CIA 
in the general operator form 

(T*(Pl,M\T t )\Aj»B\T 1 (P t ,M,T 1 ,)) , (B26) 

where is the operator of a single-nucleon current, and 
A and B are operators composed of off-shell propagators, 
on-shell projectors, and, in the cases of diagrams (D), 
(E), and (F), permutation operators. We can now insert 
a complete set of partial wave states (|B1[) before A and 
after B in (|B26[) to obtain the desired current matrix 
elements in partial wave form. 

Note that in all diagrams, as part of the states (jBljl . 
single-nucleon momentum kets \pi) that are inserted on 
the right of B are contracted to corresponding bras from 
the states inserted on the left of A. These inner products 
lead to Dirac delta functions that insure momentum con- 



servation along the nucleon lines in each diagram. In the 
following sections, we show how these delta functions are 
evaluated, and we display the final form of the various 
diagrams in partial wave form. 



1. Calculation of diagram A 



The bound state vertex vector represented on right 
(left) side of diagram A is an integral over momen- 
tum variables q, p (q' , p') and angular variables 0, 
(j), 8 ((!>', 0' , 0', 9'), and its total momentum is Pt 
{PI = B( V e 3 )P t ). 

The product of momentum conserving Dirac delta 
functions in diagram A is: 



)v(m, 0) — Pup t e ,0-Rtt ,7t ,oB(£,(q)e )v(m, 0) J 
x S 3 (B(r)e 3 )^ ie ',oZ(q')R^ >§ , fi B(((p')e 3 )v(m, 0) - R^ >e flZ(q)R $Ao B^{p)e 3 )v(m, 0)) , (B27) 

where we introduced the auxiliary four-vector v(x,y) = (x,0, 0, y) and defined sinh£(.T) = x/m. 

From the first Dirac delta function in Eq. (|B27[) we get = <l>, 0' = w — @i, q 1 = rj., where 0\, r\ are obtained by 
applying (JB9J) to B{- V e 3 )R ^-e,oB{t{q)e 3 ): 

B{- V e 3 )R^~evB{ti{q)e 3 ) = i?o, ei ,o5(£(ri)e 3 )i?o,e 2 ,o ■ (B28) 
To evaluate the second delta function, we begin by applying rule (|B9p twice, 

Z{-q)R a ^efiB{f]e 3 ) = flo^o^Si^-Ro.^.o , (B29) 
B{!i{ Sl )e 3 )Ro,o 2+ e< fi Z{q') = Ro^oB^(s 2 )e 3 )R OM , (B30) 

which defines s\, S2, J?2> $3, and #4. 

Next, we change the angular integration over <\>' and 9' to the new angles <fi and 9 (with < <fi < 2ir and < 9 < it), 
where </>' and 9' are given in terms of </> and 9 through the application of (|B17[) , 



and we have 



Using (|B9[) in the form 



#0,-^4,0^,0,0 - Rfrfi'-ip > ( B3 1) 



/ d9'sm9' dfi = / d8 sin 9 dep. (B32) 
Jo Jo Jo Jo 



B(as2)e 3 )R j fi B(ap')e 3 ) = iWo£(£(s 3 )e 3 )#o,tf 6 ,o . (B33) 
determines S3, and finally applying again (|B17I) . 

-^O,i5i+i? 3 ,0#O,i9 5 ,O = R<f>x,® 7 ,fa > (B34) 

defines the angles 4>i and $7. 

What we have achieved is that the first momentum in the second delta function has been rewritten as 

B(rie 3 )Rv,e',oZtf)Rp t9l>0 B(£tf)&)v(m, 0) = Rp,e l oZ{q)R <t>1 ^ r , B^{s 3 )e 3 )v{m, 0) , (B35) 
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where we have also used the fact that rotations leave v(m, 0) unchanged. Comparing (|B35[) to the second momentum 
in the argument of the same delta function, we can read off <f> = <j>i, 9 = #7, and p = S3. 

Before writing down the final expression for diagram A, we have to address the issue of what consequences might 
arise from the limitation of the spectator momentum in the vertex function to values below the critical value q s . The 
photon momentum absorbed by the 3N system in the lab frame is directed along the positive z-axis, which imparts 
on it a corresponding boost into the same direction. It is clear that even if the spectator momentum of the 3N system 
at rest lies below q s , it can exceed the critical value after being boosted. Whether this is the case or not depends on 
the size of the boost, the initial spectator momentum, and its angle with respect to the direction of the boost. 

We define K to be the z-component of the photon momentum in the lab frame, which is related to the invariant 
square of the photon four- momentum, —Q 2 , through K = Q^l + Q 2 /4M?. Then the boost rapidity r\ is given 
through sinh.77 = K/M t . An analysis of the effect of this boost on the spectator momentum shows that, for a given K 
and a fixed q < q s , the magnitude of the boosted spectator momentum q' is less than q s if and only if G I, where 
/ = [0, 7r] if < q < a and I =] arccos b,ir] if a < q < q s , with 



(M 2 - m 2 )y/M? + K 2 - (M 2 + m 2 )K 
2M 2 



M tv /m 2 + q 2 - y/(m 2 + q 2 )(M 2 + K 2 ) 
~q~K 



(B36) 
(B37) 



These relations are valid as long as < K < (M 2 — m 2 )/2m, which is always the case in the calculations of this work. 
Note that q' — > q s as — > arccos & when a < q < q s . We restrict the integration over to the interval /. 
Diagram A can now be written as an integral over q, 0, p' , 9, <f>, <P: 



1m 2 f qs a 2 f f + °° v' 2 f 27T - f 27r 

( ** |J!M =wl d «mL d0eine l dp 'w)i SeinS l d *l * £ 



T j m X 1 X / 2 X' 3 p 3 
TjmX 1 X 2 X 3 p 3 



VWTT^JTTv^l,^, (<*>', 0', o)v^ 2 >_ Xi (*, 0, o)v% 2 _ x , tf, 9', o) 

x Ciq'E^p'M'j'm'X'^X^ + p' 3 T l T z )C{qE{p)pMjmX 1 X 2 X 3 + p 3 TT z ) , (B38) 



where the following spinor matrix elements appear: 

F xX =u + (0,X[)S-\B(Z( q >)e 3 ))S-\R^^o)S-\R ,e>:o)S-\B( v ^ 
xS(B(£(q)& 3 ))u + (p,\i), 
FW^u+^X'^S-^BiZtf^S- 1 ^ 

x S(Ro, ,o)S(Z(q))S(R^ § )S(B(ap)e 3 ))u+(0, A 2 ) , 

F <M» =JH ™ \ uP'Hi' X') 

PkKt*** E(p') E(p) l(p< 3 + l)E(p>) - W(q')][(p 3 + l)E(p) - W(q)} 1 ■ ' 



x5- 1 ( J R 7r , 7r , )5- 1 (^, ie -, >0 )5- 1 (Z(g'))5- 1 (iio,e',o)S- 1 ( J B(r ? e 3 ))7 I/ 
x SiRo^^SiZiq^SiR^^SiR^^o^ip, X 3 ) , 



(B39) 
(B40) 



xS-\R^,a)S-\R^j, :0 )S-\Z(q>))S-\Ro,e,,o)S-\B(T ] e 3 ))YS(R ,e,o) (B41) 
x S(Z(q))S(R } ^^(i^oK 3 ^, A3) , 

F C33)^r = _m m 1 , 

p's^p^ E(p')E(p)[(p' 3 + l)E(p')-W(q')][(p 3 + l)E(p)-W( q )} [P ' 3> 
x S-^R^^^S-^R^^^S-Hziq'^S-'iRo^^S-^B^e 3 )^ 
x S(Ro,0,o)S(Z(q))S(R l§ Q )S(R ni ^ o )u p3 (p, A3) • 

c' 17— i 



(B42) 



(B43) 
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Model f(p' 2 ,p 2 ) g(p' 2 ,p 2 ) Mp' 2 ,P 2 ) F Sp (Q 2 ) F 3n (Q 2 ) 
NCI 1 1 

NCII Mp' 2 , P 2 ) Mp' 2 ,p 2 ) go( P ' 2 ,p 2 ) G Ep (Q 2 ) G E n(Q 2 ) 

NCIII fo(p' 2 ,p 2 ) fo(p' 2 ,p 2 ) 90 (p' 2 ,P 2 ) F lp (Q 2 ) F ln (Q 2 ) 

NCIV Mp' 2 ,p 2 ) 1 go(p' 2 ,P 2 ) G Ev {Q 2 ) G e „(Q 2 ) 

NCV Mp ,2 , P 2 ) 1 go(p' 2 ,p 2 ) F lp (Q 2 ) F ln (Q 2 ) 

TABLE VI: The off-shell nucleon form factors /, g, h, and the electromagnetic nucleon form factors F3 P and F 3n in ()B44[) to 
(|B46I) for the models of the nucleon current used in the numerical calculations of this work. 



Expression (|B38[) also contains the isospin-dependent functions 



4%t. = M><&) ( F lp (Q 2 ){T%\^-*-\TT z ) + F ln (Q 2 )(T% 



F^ TT =~ 9 (q' 2 ,q 2 )[F 2p (Q 2 ){T'T z 



TT Z ) + F 2n (Q 2 ){T% 



F^tt. = ( F 3p (Q 2 )(T'T z \ 



TT Z ) + F 3n (Q 2 ){T'T z 



1 - T 3 \ 


(B44) 




(B45) 


l -^\TT z )) „ 


(B46) 



where t| is the isospin projection of nucleon 3, and the four-momenta q 3 and q 3 are defined as 

% = (W(q') - E(p'),0,0,p') , q 3 = (W(q) - E(p), 0,0, p) . 



(B47) 



We have constructed five different models of the electromagnetic nucleon current, which correspond to the different 
choices of the off-shell nucleon form factors given in Table IVT1 
The isospin matrix elements needed in (|B44[) to (|B46[) are 

(0T z |i^i|07;> = i, { 1 % \ 1 -^\1%) = 1±1%, <lT z ' 1±r3 



1 



1±T. 



F { f^, v , in (IB4T1) and F ( , 3 ^ l/r . in (lB42l) contain factors of the form 



(B48) 



Eft) (P3 + 1)E[ p)-w( q ) ( and the analogous 
expression with primed variables). These factors originate from the p-spin decomposition (| 1 5[) of off-shell nucleon 
propagators. We introduce 



m 



i 



E(p) (p+l)E(p)-W(q) 



(B49) 



as the p-spin propagator of nucleon 3. 

The positive p-spin propagators g + (q,p) and g + {q',p') are finite since 2E(ft') — W{q) > 3m— Mf > 0. On the other 
hand, the negative p-spin propagator develops a singularity as q approaches q s . In the vicinity of q s , 



Similarly, 



9 (q,p) 



g-(q',P') 



1 



1 



W{q) y/%fo(«.-«) 



W(q') y/2Kq(b-cosG) 



(B50) 



(B51) 



as O — > arccos6 and a < q < q s . These propagators are only weakly singular, and they are still multiplied by the 
vertex function which goes to zero at the singularity. The total result is finite. 

The integrations over $ and <j> in Eq. (|B38|) can be done analytically. The only ^-dependence in the integrand of 
(IB381) resides in the rotation matrix (i?<z> : o,o) M „ and in a factor e -m'-M)4> from 2?j^, (<£', 0' , fyV^m-Xi & °> °)- 
The product is easily integrated: 



( * dQe-^'-^Rzw = (2 - \M' - M\)ir 
Jo 



( 1 - \M' - M\ \ 

\M'-M\ i{M'-M) 

-i(M'-M) \M'-M\ 

\ 1 - \M' - M\ J 



(B52) 
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Next, the integrand of (|B38|) depends on <j> through </>', 9', <fi, 9, which appear in the factors V ( £\,_ y (4>',9',0), 
D^Dt , (d>, 9, 0), -Fw , , F^f 3 }," . . F < '-f 3 } VT . and G'fl" , . We can write those spinor matrix elements as 

m,A 2 -A 3 v^' ' n A 2 A 2 ' P3A3P3A3' P3A3P3A3 P3A3P3A3 * 

f kI= E E ^f(^^o)p<^ 

S6,=±i/Jt^=± (B53) 
x [u T ^0,X 2 )u+(p,X 2 )}u<(O i e 2 )S-\Z(q'))S~\^e', )S-\B( V e 3 ))^ 

ViM yP> ? 3 ?3=±l/2^T3 = ± 

(B54) 

x^o^S-^Z^S-^iW.o)^ 1 ^^ As)], 

(B55) 

< ; ZU = E E ^W^K 1 ^ 

W ? 3 C3=±1/2t 3 T3^3 = ± 

xu^(0,QS-\Z(q'))S-\R ,e',o)S-\B( V e 3 ))YS(Ro,e,o)S(Z^ X 3 ){u^(0,X 3 )u T *(p,X 3 )}. 

(B56) 

Hence the integration variable <fr appears only in the product 

V { £] K ^\e\0)V%l 2 ^ (B57) 
which can be written as 

-/ 

s'=-j' s=-j /%/3 3 fts/3 3 =±l/2 

x d££ + tf 3 )d£S (WSfc + ^S-l ^ ■ (B58) 



The (f> dependence is now isolated in the factor e li - s +l 3 2+03+ s +l 3 2+M<t' ; and the analytical integration over has 
become straightforward. 



2. Calculation of the sum of diagrams B and C 



Diagrams B and C are closely related by symmetry. We start with diagram C by evaluating its momentum- 
conserving Dirac delta functions, 



6 3 (B(r,e 3 )R&' ,0' fiR-K,TT : oB{^{q')e )v(m, 0) — R$,e,oRn,ir,oB(£(q)e )v(m,0)) 

x 8 A ( K B{^)R^ l ,e'fiZ{q')R^ s , >Q R^^ Q v{W{q') - E(p'),p') - R^,efiZ{q)R u>Q R v ^ Q v{W{q) - p°,p)) • (B59) 

The first delta function in (|B59[) is the same as the one in (|B27|) of diagram A. Therefore = 0' = tt — &x, 
q' = r\ , and the interval / for the 0- integration is also the same as in the case of diagram A. 

We use the variables #3, s 2 defined in Eqs. (|B29j) . (|B30|) . and <[>, 9, cf>' , 9' which are related through (|B31|) . Now 
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we introduce S3 and $5 (with < $5 < 7r), 



S3 



\ 



[p' sin d) 2 + f [WW) - E{&] ) — - p> cos 6 + ) , (B60) 
\ m ml 



s 3 cos i? 5 = - [WW) - £(p')l — + P' cos 0-¥ ?- . (B61) 

m to 

Next we define </>i, i9g , </> 2 , S4 by 

^0,191+133,^-^0,195,0 = ^l,l? 6 ,<^2 ! (B62) 



s 4 = W(g) - [WW) - £(£')] + S ' + p cos 0— . (B63) 

TO TO 

Evaluating the second Dirac delta function in Eq. (|B59p in the way it was shown for diagram A, we obtain p a = S4, 

p = s 3 , 4> = 4>i, 6 = j? 6 - 

Diagram C can now be written as an integral over q, 0, p' , 9, <f>, <P: 

3to 2 f qs a 2 f f +co fi' 2 f 27r - f 27r 

^■swi d "W)l deei " e l *mih*'* l k di L d \, , £ 



T'j'm'A'jA^Agp^ 
TjmXi A 2 A3P2P3 



v/2^y2j+T^,_ A , (<*>', 63', 0)D^: Ai (*, ©, 0)Z>^ A , _ A , (0', 0', 0) 

x 0)C7(< Z 'S(p')p'M'j , ™'A 1 A^3 + p' 3 rT z )C(qp pMjmX 1 X 2 X 3 p 2P3 TT z ) 

x ( W, (4^£a 2 + ^S>) <?*U , (B64) 



where Fy 1 ? is defined in Eq. (|B39|) . The remaining spinor matrix elements are 



A' x Ai 

F [ 22 > = P° + P 2E }P} m n+(p', X^S-^Rr, g , JS-^Ztf^S-^Ro e , 0) 



x 5- 1 (B(7 7 e 3 )) 7 ^(i?o,e,o)5(Z(g))5(i?^ K 2 (p, A 2 ) , 



x S-\B(r ) e i ))ia VT S(R o , ,o)S(Z(q))S(R^ lo )u'''(p, A 2 ) , 
,(33) _ m , , , ( to 



xS- 1 (fl T ,,,o)S- 1 (fl|, f , )S- 1 (2( g , ))S- 1 (^,e',o)S- 1 (5(')e 3 ))S(fio,e,o) 
x S(Z(q))S(R $AQ )S(^,o)u^(p,X 3 ) , 



(B67) 
(B68) 



and the isospin dependent factors are 



F^ TTz = f(m 2 ,q 2 2 ) (F lp {Q 2 )(T'T z \ l -±^\TT z ) +Fln (Q 2 ){T'T z \ l -^\TT z )j , (B69) 

4%% = 9(m 2 ,q 2 2 ) (V 2p (Q 2 )(T'T 2 |i^i|TT 2 ) + F 2n {Q 2 )(T'T z \ l —^\TT z )^ , (B70) 

with 

ga = (p°,0,0,p). (B71) 

The matrix elements involving the isospin z-projection of nucleon 2, r|, are 

(or,|i^i|or 2 ) = i, ( ir 2 |i^i|ir 2 ) = i±^, ( it 2 |I^|ot 2 ) = T ^r 2 , (or 2 |i±i|ir 2 ) = T ^r 2 . 

(B72) 
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Just as it was shown for (|B38[) . the integrations over ^ and <fi in (|B64|) can again be done analytically. 

We turn our attention now to the integration over 9. This is complicated due to the moving singularity contained 
in the factor (p° + p2E(p))/(E 2 (p) — (p ) 2 ) in (|B65j) and (|B66|) . which has its origin in the off-shell propagator of 
nucleon 2. In this work, we calculate diagram C separately from B, but eliminate the imaginary part of the propagator 
singularity, which cancels with the corresponding one of diagram B, and keep only the principal value integral. The 
energy denominator can be written as 

E 2 (p) - (p ) 2 = AE(p') - Bp'cos9-C, (B73) 

where 



A = 2W(q')-2W(q) ^ m2 + S2 , B = 2W(q) S -l, C= A \ B2 . (B74) 

m m 4 

A, B, and C do not depend on p' nor on 0, and B > \A\. The location of the singularity depends on the relations 
between A, B, and p' , and we have to distinguish between several cases. We define 



and, for p' ^ 0, 



BC±\AWC 2 + m 2 {B 2 -A 2 ) 
P± = ^ (B75) 



We find that the equation E 2 (p) — (p ) 2 = in the unknown 9 has at least one solution if and only if one of the 
following conditions is satisfied: 

A > 0, p' > p- , (B77) 
A < 0, B < yfA 2 + Am\A\, p' > -p+ , (B78) 



A < 0, B = y/A 2 +4m\A\, p > , (B79) 



A < 0, B > y/A 2 +4m\A\, p'>p+, (B80) 
p=0. (B81) 

The respective minimum values p_, ~p+, and p + of p' in (|B77|) . (|B78[) , and (|B80[) are all positive. In any of the cases 
(|B77|) - (|B80j) . the equation E 2 (p) — (p ) 2 — has only one solution, namely 9 = arccos D(p'). In the case of (|B81j) . 
the set of solutions of E 2 (p) — (p ) 2 = is the interval [0, ir]. 

Now that the location of the singularity has been determined, one can divide the integration region into subintervals 
according to the four cases (|B77[) - (|B80I) and apply standard numerical techniques (we used a subtraction method) 
to perform the principal value integration. 

Diagram B does not need to be calculated once diagram C is known, since it can be obtained from the latter using 
time reversal and parity transformation properties. It can be shown that the components of diagram B are given by 

(M'\J° B \M) = {B(^))\{M'\4\M) - (B( V e 3 ))%(M'\ J*\M) , (B82) 

(M'\ Ji\M) = (MVcW . ( B83 ) 

(M'\J%\M) = (M'\ J 2 C \M) , (B84) 

(M'\J%\M) = -(B(r,e 3 )) 3 3 {M'\4\M) + (B( V e 3 )) 3 (M'\J c \M) . (B85) 

3. Calculation of diagram D 

The product of momentum conserving Dirac delta functions in diagram D is 

S 3 (B(r]e 3 )^,0,, o Z(q')R^ iS , !O B^(p')e 3 )v(m, 0) - R^,0 tO R^^ o B^(q)e 3 )v(m, 0)) 

x S 3 (B(i 1 e 3 )R <P ,, e , fi R^ 7rfi B{Z {q')e 3 )v (m, 0) - R^efiZ(q)R^ § B^(p)e 3 )v(m, 0)) . (B86) 
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From the first Dirac delta function in (|B86[) we obtain <f>' = (ps, 8' = and p' — s 2 , where (fa, #6j S2 are defined by 
applying (JB9J) and ([BIT)) : 

Z(-q')i?o,-e',oS(-7]e 3 ) = i?o^ 1 ,oS(e(si)e 3 )i?o,^ 2 ,o , (B87) 

Ro,d 2 ,oR&-<P'+Tr,TT-0,O = Ripi,$ 3 ,ipi : (B88) 

B(Z( Sl )e 3 )R ^oB(Z(q)e 3 ) = i?o,,? 4 ,o5(£( S2 )e 3 )i?o,,? 5 ,o , (B89) 

•Ro,i?i,yi-Ro,i?4,0 = Rip 3 ,-8(i,ip4 ■ (B90) 

The second delta function leads to <p = fa, 9 = 8 6 , and p = r 2 , where — again by using (|B9|) and (|B17|1 — fa, 9q, r 2 
are determined through 

Z(-q)Ro t _ 0fl B(r]e 3 ) = i?o, ei ,oS(^(r 1 )e 3 )i? , e2 ,o , (B91) 
Ro,9 2 fiR<p'-<p+Tr,Tr-0',o = R<t>t,e 3 ,<j>2 > (B92) 
B(ari)e 3 )Ro,e 3 , B(aq')e 3 ) = i?o,0 4 ,o£(£(r 2 )e 3 )i?o, e5 ,o , (B93) 
Ro,e 1 ,oR<j>i,e i ,o = R<p 3 ,e 6 ,<f>4. ■ (B94) 

We can write the current given by diagram D in the form 

3m 2 f qs a 2 I" 13 a' 2 r r f 27r f 2 * 

www—wl *mi «mh de * me l dew i *7. * 

2 ^ v^^>/Srn^.^(#'.^.o)i>^^(#,e,o)i^!» s . Al (i',#'.o) 

TjmXx A 2 A 3 p 3 

x ^i-AaC^^^^C^^^Op'^'j'^'AxA^Ag + p' 3 T'T z )C(qE(p)pMjmX 1 X 2 \ 3 + p 3 TT z ) 

F (12) p (21) / R \" p (33> p (2) p (33)^r p (3) r (33> \ mq ,^ 

X AJA 2 ijii ^£+£i, ,oJ „ (/VtT^A^As + T'TT Z p' 3 X' 3 p 3 X 3 It + b T'TT^ p' 3 X' 3 p 3 X 3 ) » (B95) 

with the spinor matrix elements 

F^ ' =w+(0,Ai)5- 1 (B(e(gOe 3 ))^" 1 (^,o)^ 1 (i?o,e',o)5- 1 ( J B( 7? e 3 ))5( J R*_^ >e>0 ) 



x S(Z(q))S(R^ § fi )S(B(aP)e 3 ))u+(0, A 2 ) , 
Fjg = fi+(0, ^S-^BiZffl&yS-HRw, )5- 1 (Z(g'))5- 1 (i?o,e',o)^ 1 (S(r/e 3 )) 
x 5(i? <f _ s ,,, e ,o)5(^^,o)^(S(e(g)e 3 ))u + (0, Ax) , 

1,(33)* 



TO „',,,, TO 



(B96) 
(B97) 



x S- 1 (i^, w>0 )5- 1 (i?^ )§ , )0 )5- 1 (Z(gO)5- 1 (iio,e%o)5- 1 (B(r ? e 3 ))5(i? £ _^ i0)0 ) (B98) 
x 7 t/ g(i?^ !O]O )g(J?o.0.o)g(^(q))g(^ ^n)^(^. ff .o)n p3 (p, A 3 ) , 

F (33>r m / , , TO -gWPsfo' X'o) 

r P ' 3 x' 3P3 x 3 E(p') 9 [q,P) E(p) 9 [q,p> yP> A z) 

x S-\R^^o)S~\R ilt§ , fi )S-\Z(q'))S-\Ro^,o)S-\B( V e 3 ))S(R^ Q Q ) (B99) 
x ia^sm.-* fifi )S{R .e.o)S{Z{q))S{R^ i5 n )S{R lt ^o)ui' 3 (p, As) , 

G ^ - (^)) 2 

x5- 1 (K ! ,o)r 1 (^,|,„)S- 1 (Z( 9 '))5- 1 (flo, ff ,o)r 1 (B(» / e 3 )) (B1 " 0) 
x S(R*-& Q0 )i v S(R ± - ±Lq Q )S(R Q .0 fi )S(Z(q))S(R^ g fi )S(RT Ty7Tfi )u' T (p, A 3 ) , 
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and the isospin-dependent factors 

n(l) _ tf„'2 „2 



F&TT. = f(q' 3 2 ,qi)(F lp (Q 2 )(T'T z \P% 
F(% % = ~g($, q 2 3 )(F 2p (Q 2 ){T'T z \P% 

4%T. = kq'i,ql)(F 3p (Q 2 )(T%\Pf°^-±\TT z ) + F 3n (Q*)(T'T z \P£ 



1 


+ r! 




2 


1 


+ r| 




2 


1 


+ r! 


2 



\TT Z ) + F ln (Q 2 )(T'T z \P± 
TT Z ) + F 2n {Q 2 ){T'T z \P{ 



1 






2 


1 


- T 3 




2 


1 


- T 3 

r 3 


2 



\T%)), 

TT,)), 
\TT Z )). 



(B101) 
(B102) 
(B103) 



The four-momenta g 3 and q 3 are defined in (|B47|) . where now q' is an independent variable, and we have p' = s 2 and 
p = r 2 . P|| i s the operator that interchanges the isospin states of nucleons 1 and 2, with the matrix elements 



(it z \p% 



1 + T 3 

3 , or : 



4 2 

4 T 6 2 



(ir 2 |p; s 2 ° 



l±r|, _ 1 5^ 
2 17 4^6 



1 + -r 3 

3 'IT, I 



2 ' ' 4 6 2 

In (|B95[) , we change now integration variables from <P' and <P to <fi and (/J, with 

(j> — (£>' <P -\- <P' 

(f>= — n — , ¥>= — o — > 



4^6 



(B104) 



(B105) 



which implies 



2tt 



2tt 



d<P' / d<P = 2 
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•.(V2) a , n^m* 



(B106) 



The only 93-dependence in the integrand comes from V^, ; nl _ x (^',0 / ,O)2?^* Ai (<^,6),O), which yields a factor 
e -i{M -M)<^p and can be integrated analytically. The three different cases are: 



if M' = M 



r -\<j>\+2TT 

J\*\ 



if M' = -M = 1/2 



-I0I+27T 



( -2sin|0| 



X 



/ 2tt — 2j0| \ 

-2sin|</>| 

—2 sin|^| 

\ 2<n-2\<j>\) 





-|0| + 7T - i sin(2|0|) »[-|^| + tt + i sin(2|0|)] 

-i[-|0| + tt + i sin(2|^|)] -101 + tt - i sin(2|0|) 





(B107) 



if M' = -M = -1/2 



-\4>\+2n 



d^e-W-^IL, ,0,0 = 



/-2sin|0| 

-101 + 7T- isin(2|0|) -i[-|0|+7T+isin(2|0|)] 

z[-|0|+7r+isin(2|0|)] -|0|+7r-isin(2|0|) 

V 0" 



\ 




-2sin|0| / 

(B108) 



\ 



-2sin|0| J 

(B109) 



As in the case of diagram A, there is no problem with propagator singularities and the numerical integrations can be 
performed in a straighforward manner. 
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4. Calculation of the sum of diagrams E and F 



In close analogy to diagrams B and C, also diagrams E and F are related by symmetry. We start with the calculation 
of diagram F and relate it later to diagram E. Here, the product of momentum conserving delta functions is 



x S i ^R4>,0, o Z(q)R^ § o v(p°,p) - B(r}e 3 )R$, l e, fi R Kl „ fi B(t(q')e 3 )v(m,Q) + q 



(B110) 



The first delta functions is the same as the one in (|B110|) . hence fa = if3, 9' = 'Qq and p' = S2- 
We define r\ and 9\ as 



\ 



sin 2 & + I —E(q') - ^S±Hq' cos 0' - K] , 



ri cos 9 l = — E(q') - ^ M t +K \ ' cos Q> - K , (0 < 0i < tt) 
M t M t ~~ ~~ 



(BUI) 
(B112) 



and introduce fa, fa, $2 by applying (|B17[) , 



-Ro,-e,-<J>+<J>'+7r Rofii,o — R<t>!, 



(B113) 



From the second Dirac delta function in Eq. (jBHOj) we obtain 



P W(q) \ W > M t q M t V * J W(q) 



v\ cos Qi 



(B114) 



as well as cj> = fa, 9 — 8 3 , and p = r 2 , where 



r-2 



t\ sin 2 6 2 



I W J Aft y M t V V i 



ri cos 02 



1/2 



(B115) 



r 2 cos 9 3 



y/W 2 (q)+r 
W(q) 



r lC os6 2 , (0<6» 3 <7t). 



(B116) 



The current matrix elements of diagram F are given by 



(M'\J$\M) = I dq q 



2tt r 2ir 



T^TS r dq^— I dB sine I d&l d^T'dq'-^-^l dO' sin 0' 
(2tt) 8 J q q E(q)J J J J q E(q')J 



T'j'm' A^ X' 2 
Tj'mAi A 2 A 3 p 2 p3 

2 — A3 v 



VWTiVWTTv^%,_ K (<*>', 0', o)v^\ (*, 0, o)v^] x ,_ x , e>, o) 



X 



^L-LfA 0, 0)C(^(pO^M'jWAiA' 2 Aj + p , 3 T'T z )C(qp»pMjm\ 1 \2\ 3P 2P3TT z ) 



X Ai ^±+£1,0,0,/ „ ^T'TT/A' lP2 A 2 +- t T'T% ± '\< 1 p 2 \s ( h-) f p < a \< 3 



3P3A3 ' 



(B117) 
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where Fy Ai is defined in (|B97|) , and the remaining spinor matrix elements are 



x S&s-p^SiZiqVSiRM K 2 (p, A 2 ) , 



(B119) 



(B120) 



x S(Re-* i0tO )S(Z(q))S(Ri ilo )v l '» (p, A 2 ) , 

x S- 1 (J2 w ,» l0 )5- 1 (^ (fl -,o)5- 1 (^(?'))5- 1 (i2o,©',o)5- 1 (B(^ 3 )) 
x S(E*_^, e ,o)5(Z(g))5(^ e - )5'( J R 7r>7r ,o)u <T3 (p, A 3 ) . 
The isospin dependent factors are given by 

4^ = /(m 2 , ? 2 2 ) ^i P (g 2 )(T'r 2 |pf 2 °i±^i|rr 2 ) + ^ ln (Q2 )(T% | P i|ol^i| rrz ^ ; (B121) 

4'U = K^ 2 , <?2) (F 2p (Q 2 )<T'T 2 |Pr 2 °i^i|rr 2 ) + F 2n (Q 2 )(T'T z |Pi 2 °l^i|TT z )^ , (B122) 

where q 2 is defined in (|B71|) . its component p° now being determined in (|B114|) . and p = r 2 . The matrix elements 
involving the isospin projections of nucleon 2 are 

(0T z \P™^-\0%) = \±\% {IT Z \P™ 1 -±^\1T Z ) = -i ± ir 2 

(ir.|ifr^|QT.) = ± m \p^-^l\ir z ) = T . (B123) 

It is again useful to change integration variables in (|B117j) from <&' and <P to and according to (|B105p and (|B106|) . 
The integration over ip can then be carried out analytically, in analogy to the integration over ip in diagram D. 

As in the case of diagrams B and C, we encounter propagator singularities, whose location needs to be determined. 
For this purpose, the denominator of the factor (p° + p 2 i?(p))/ (E 2 (p) — {p ) 2 ) which appears in (|B118p and (|B119[) . 
is written as 

E 2 {P) - (p ) 2 = Q 2 - 2Q(sinh(C + v)E(q') - cosh(C + vW cos O') , (B124) 

with 



K 2 - O 2 

sinh C = - A / — TvT~ ■ (B125) 



We find that E 2 (p) - (p ) 2 = if and only if a < q' < q s and & = 0' s , where 



<^-- n )E[Q) -sinh(C + 



(B126) 



sinh(C + r?)^') - § 



B(g') = ^ , ,/ \ 2 , ©I = arccosB(g') . (B127) 

q' cosh(C + rj) 

With this information it is now possible to carry out the principal value integration over 0' . The remaining integrations 
are straightforward. 

Finally, diagram E is obtained from diagram F in complete analogy to the way diagram B is calculated from diagram 
C, namely 

(M'\J° E \M) = (B(r,e 3 ))° (M'\J° F \M) - (B( V e 3 ))%(M'\J F \M) , (B128) 
{M'\J^\M) = (M'\J F \M) , (B129) 
(M'\J 2 \M) = (M'\J 2 \M) , (B130) 
(M'\J%\M) = -{B^))\(M'\ J%\M) + (B( V e 3 ))\{M'\J F \M) . (B131) 
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